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ABSTRACT 

We present 1 6 GHz (1.9 cm) deep radio continuum observations made with the Arcminute 
Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine 
these new data with archival information from an extensive literature search to examine spec- 
tral energy distributions (SEDs) for each source and calculate both the radio and sub-mm 
spectral indices in two different scenarios: (1) fixing the dust temperature (7^) according to 
evolutionary class; (2) allowing Tj to vary. We use the results of this analysis to place con- 
straints on the physical mechanisms responsible for the radio emission. From AMI data alone, 
as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of 
the objects in this sample have spectral indices consistent with free-free emission. We find an 
average spectral index in both T^ scenarios, consistent with free-free emission. We examine 
correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow 
force and find that these data are consistent with the strong correlation with envelope mass 
seen in lower luminosity samples. We examine the errors associated with determining the ra- 
dio luminosity and find that the dominant source of error is the uncertainty on the opacity 
index, j3 . We examine the SEDs for variability in these young objects, and find evidence for 
possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1. 

Key words: radiation mechanisms: general — ISM: general — ISM: clouds — stars: forma- 
tion 



1 INTRODUCTION protostellar evolution (Andre et al. 1993) in which the protostar 

has yet to accrete most of its mass from the surrounding envelope, 
Young stellar obiects (YSOs) are divided into a number of evolu- j »i r ^i . i i ■ . ■ i ■ .i_ .i i 

'^ J V / gjj(j thereiore the central object is less massive than the envelope 

tionarv classes. Class protostars represent the youngest phase of ,,, ,. n r,,, , ^, t t /t i .r.r,-,x 

-' ^ ^ J b f (Menv > M*). The subsequent Class I phase (Lada 1987) occurs 

when a majority of the remainder of the envelope has accreted onto 
the protostar or its circumstellar disc (Menv < ^*)- The Class II 

* We request that any reference to this paper cites "AMI Consortium: ^^^^ (Lada 1987), also known as the Classical T Tauri stage, oc- 

Ainsworth et al. 2012". 
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curs when infall is almost complete and the central object is em- 
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bedded in an optically thick disc and enters the pre-main-sequence 
(PMS) stage of evolution. Embedded YSOs are often found to have 
partially ionised and coUimated outflows/jets driven by the accre- 
tion from their envelopes, and are detectable from molecular gas 
tracers such as CO and by the shocks along the length of the out- 
flows. The outflows driven by Class I protostars are less powerful 
than those of Class due to a decline in the accretion rate (Bon- 
temps et al. 1996). 

The target sample in this work is comprised of Class and 
Class I protostars. Since the mass ratio Menv > -M* is not directly 
observable, Andre et al. (1993, 2000) define the following observa- 
tional properties required to identify an object as Class 0. The first 
is the detection of a compact centimetre radio continuum source, 
a collimated CO outflow, or an internal heating source as indirect 
evidence for a central YSO. These criteria distinguish the central 
object from a starless prestellar core. The second requirement is 
centrally peaked but extended submillimeter continuum emission 
tracing the presence of a spheroidal circumstellar dust envelope (as 
opposed to just a disc). The third is a high ratio of submillimeter to 
bolometric luminosity, corresponding to an envelope mass greater 
than the stellar mass. Class spectral energy distributions (SEDs) 
are typically well characterized by single temperature blackbodies 
with 15 ^ Td ^ 30 K and where Tj is the dust temperature. The sec- 
ond and third criteria differentiate Class sources from the more 
evolved Class I and II objects. For further distinction, bolometric 
temperature is also used to characterise evolutionary class, where 
Tboi < 70K for Class objects and 70 < Tboi < 650 K for Class I 
(Chen et al. 1995). 

The radio counterpart at centimetre wavelengths for these 
low-mass young stars is detected almost certainly via thermal 
bremsstrahlung radiation. This "free-free"emission is typically ob- 
served to be compact and elongated along the outflow axis, lead- 
ing to the designation of these sources as '"thermal radio jets"'(e.g. 
Rodriguez 1995). Free-free emission at these wavelengths typically 
has a flat or positive power-law spectral index a, where the flux 
density Sy °= v" at frequency v. It has been shown that a = 0.6 
for a standard conical jet, and a < 0.6 for an unresolved, partially 
opaque flow where the cross-sectional area grows more slowly than 
length (Reynolds 1986; Anglada et al. 1998). However, in a num- 
ber of cases non-thermal emission is also seen (e.g. Ray et al. 1997; 
Carkner et al. 1997; Feigelson et al. 1998; Carrasco-Gonzalez et al. 
2010a). There are some objects that have yielded spectral indices 
too negative to be explained by free-free emission alone, such as 
the Serpens MMS 1 triple radio source (Rodriguez et al. 1989b; 
Curiel et al. 1993; AMI Consortium: Scaife et al. 2012b), suggest- 
ing that non-thermal processes must contribute. The radio emission 
can also be attributed, at least at very short wavelengths, to the discs 
around these young stars (Rodriguez et al. 2008b). 

Here we present 16 GHz observations of a sample of classic 
low-mass young stars driving outflows. These sources are selected 
to match the target sample for the e-MERLIN (extended Multi- 
Element Radio Linked Interferometer Network) legacy project at 
5 GHz on the morphology and time evolution of thermal jets as- 
sociated with low-mass young stars (RI. Rodriguez). The high spa- 
tial resolution observations provided by the e-MERLIN programme 
will resolve out the larger scale emission from these objects and 
only detect the small-scale localized emission from the central 
source and shocks along the length of their collimated outflows. In 
contrast, the data presented here will not resolve the separate com- 
ponents of the radio emission but will measure the integrated radio 
emission from these objects. This total radio emission at 16 GHz 
from YSOs has been demonstrated to show clearly defined trends 



with the other global properties of these systems, such as bolomet- 
ric luminosity, envelope mass, and outflow force (AMI Consortium: 
Scaifeetal. 2011a,b, 2012b). If such free-free emission is to be rec- 
onciled with ionization mechanisms arising from the impact of the 
molecular outflow or protostellar jet/wind on its surroundings then 
it is vital to measure it in its entirety to determine quantitatively its 
physical correspondence to other global properties of the protostel- 
lar system, in much the same way that radio emission from massive 
young stellar objects can be used to infer the total ionizing flux of 
the central object and hence its spectral type. Such indirect mea- 
sures of protostellar activity in the early embedded phase of YSO 
lifetimes may provide a window into the evolution of these objects, 
which is otherwise obscured. 

In this paper we investigate the properties of the integrated ra- 
dio emission from this sample of low-mass YSOs. In general, these 
new data confirm trends found in previous works and improve upon 
the statistics. We identify outliers which may provide useful tar- 
gets for further study at high spatial resolution. We combine our 
results for the flux densities with those found in an extensive liter- 
ature search to calculate the spectral indices at radio wavelengths, 
and which are overwhelmingly consistent with free-free emission 
as the mechanism for the radio emission. 

The organization of this paper is as follows. In §2 we describe 
the sample of targets to be observed. In §3 we describe the AMI 
telescope, the observations, and the data reduction process. In §4 
we comment on the results of the observations and §5 contains de- 
tailed notes on the individual fields. In §6 we discuss the detections 
and non-detections, the expected contamination from extragalactic 
sources, the spectral energy distributions, the radio spectral indices, 
the derivation of physical parameters, correlations between the ra- 
dio luminosity and other global properties, and finally, evidence for 
variability. We make concluding remarks in §7. 



2 THE SAMPLE 

This study targets eleven YSOs driving known outflows. The coor- 
dinates of each of these objects are listed in Columns [3] and [4] 
of Table 1 along with their identifier. We follow the classification 
scheme described in Hatchell et al. (2007b) for Class or Class I 
protostars, based on three evolutionary indicators, which for a 
Class I are Tboi > 70 K, Lbol/ismm > 3000, and Fs.e/Fgso > 0.003, 
where F3 g is the flux at 3.6 cm and and Ff^^Q is the flux at 850 flm. 
There are eight Class and three Class I protostars in the target 
sample, listed in Table 1. A number of additional sources were de- 
tected within the AMI fields. We identify these additional sources 
assuming a generous threshold of lOarcsec of a known protostel- 
lar source to allow for low signal to noise, however the maximum 
offset is found to be < 4arcsec. Additional sources are listed in 
Table 2. When combined with the original target list, these addi- 
tional sources will be refered to as the "extended sample". For the 
purposes of this work, all other detected sources are considered to 
be extra-galactic and this is discussed further in §6. Physical prop- 
erties of the sources including parent cloud association, distance, 
bolometric temperature, bolometric luminosity, envelope mass, and 
outflow force are listed in Table 3. We note that the bolometric lu- 
minosities of the target sample span an order of magnitude, and the 
range of envelope masses places the sample in the intermediate to 
high end of the range for pre-main-sequence stars. 

A number of the physical properties listed in Table 3 depend 
on the distance to the source. For example, many studies assume a 
distance to Perseus of 320 pc (e.g. Hatchell et al. 2007b), but here 
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Table 4. AMI frequency channels and primary calibrator flux densities mea- 
sured in Jy. 



Channel No. 
Freq. [GHz] 



4 
14.62 



5 
15.37 



6 
16.12 



7 8 

16.87 17.62 



3C48 
3C147 



1.75 
2.62 



1.66 

2.50 



1.58 
2.40 



1.50 
2.30 



1.43 
2.20 



we assume a distance of d — 250 pc in order to remain consistent 
with the Spitzer "Cores to Disks" (c2d) catalogue of deeply embed- 
ded protostars (Dunham et al. 2008). Values in the literature also 
differ for the distance to the Taurus and Orion molecular clouds, 
and the values we adopt, along with their references, are listed in 
Table 3. Physical quantities from the literature have been adjusted 
for distance where necessary. 



3 OBSERVATIONS 

The data in this work were taken with the Arcminute Microkelvin 
Imager (AMI) Large Array at the Mullard Radio Astronomy Obser- 
vatory near Cambridge, UK (AMI Consortium: Zwart et al. 2008). 
AMI is comprised of eight 13 m dishes, operating between 13.5 and 
17.9 GHz with eight 0.75 GHz bandwidth channels. The channels 
1-3 were not used due to a poorer response in this frequency range 
and interference due to geostationary satellites. The FWHM of the 
primary beam of the AMI is ~ 6 arcmin at 16 GHz. 

AMI data reduction is performed using the local software tool 
REDUCE. This applies both automatic and manual flags for interfer- 
ence, shadowing and hardware errors. REDUCE Fourier transforms 
the lag correlator data to synthesize frequency channels and per- 
forms phase and amplitude calibration before output to disc in uv 
FITS format. Flux (primary) calibration is performed using 3C48 
for all fields except for LI25 1 where 3C147 was used, and assumed 
flux densities can be found in Table 4. Phase (secondary) calibra- 
tion is carried out using the bright point sources listed in Table 1. 
Absolute calibration is typically accurate to 5 per cent but we note 
that an increased absolute calibration error of 10 per cent is applied 
to flux densities from channel 8 in order to reflect the poorer phase 
stability of this channel relative to the others. 

Reduced data were imaged in CAS A'. Multiple data sets 
were concatenated, additional baseline flagging was carried out, 
and the CLEAN task was used for deconvolution to produce the 
combined frequency image as well as the separate spectral chan- 
nel images. Primary beam correction and fitting Gaussian models 
to radio sources in order to extract flux densities were performed in 
AIPS^. All errors quoted are 1 o. 



4 RESULTS 

Maps were made using naturally-weighted visibilities to ensure op- 
timal signal-to-noise levels, except in the case of HH 1-2 in which 
uniform weighting was used to improve resolution. As HH 1-2 
is the lowest-declination source, the beam would otherwise have 
been very extended. The dimensions of the synthesized beam for 
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each source are listed in Table 1. Source detection for these data 
was performed in the un-primary-beam-corrected maps, where we 
identified all objects within the FWHM of the AMI primary beam 
with a peak flux density > 5c!'r[Y,s, where tTrni.s is the rms noise de- 
termined from the maps. CTfms for each target field is listed in Ta- 
ble 1. Integrated flux densities for detected sources were then ex- 
tracted from the primary-beam-corrected maps using IMFIT and are 
listed in Table 5. The errors on the flux densities are calculated 
as (7 = J {O.OSSv)^ + (Tj? + CTf^s, where (7^ is the fitting error re- 
turned from IMFIT and 0.055v is a conservative 5 per cent absolute 
calibration error (except for channel 8 where we use a 10 per cent 
absolute calibration error, 0.15v). Including a^ provides an over- 
estimation to the total error. 

Extracted flux densities are listed in Table 5, and column [8] 
of Table 5 lists the spectral indices ojami calculated over the AMI 
channels. The uncertainties of these values are large due to the short 
frequency coverage of AMI. Spectral indices are discussed further 
in §6. In two cases, HH 26 IR and NGC 2264, the flux density 
was extracted only from the 16 GHz combined-channel image. For 
these objects, the flux density in the combined-channel image is 
< 5y^a, where n is the number of combined channels, indicating 
that the sources would be detected at < 5cT in the individual channel 
maps. Combined-channel maps for each source may be found in 
Appendix A with the exception of L1448 IRS 3 which is shown in 
Figure 1(a) for illustration purposes. 



5 NOTES ON INDIVIDUAL FIELDS 

L1448. L1448 IRS 3 is made up of the L1448 North group of 
Class protostars, NA, NB, and NW (Curiel et al. 1990; Terebey 

6 Padgett 1997). AMI does not resolve the separate sources, so the 
data at 1.9 cm (see Figure 1(a)) represents the integrated radio emis- 
sion from the three protostars and appears point-like. The spectral 
energy distribution of this source is shown in Figure 1(b), and in- 
cludes the observed radio spectra over the AMI frequency channels 
3-8 combined with flux densities from the literature. Spectral en- 
ergy distributions and spectral indices are discussed in §6. 

LI 448 NA and NB are in a 7" separation protobinary system 
(Curiel et al. 1990) and L1448 NW lies approximately 20" north- 
west of L1448 NA (Terebey & Padgett 1997) in the Perseus molec- 
ular cloud. The outflows in this region have been discussed exten- 
sively by Wolf-Chase et al. (2000). Only redshifted gas associated 
with the outflow of L1448 NA is detected in their CO maps, with a 
PA (measured from north through east) of ~ 150° and a total length 
of at least 0.7 pc (at their assumed distance of 300 pc). The outflow 
associated with L1448 NB has a PA r^ 129° and 6' to the northwest 
along the outflow axis lie the blueshifted optical emission knots 
of HH 196 (Bally et al. 1997), the directions of these outflows are 
illustrated in Figure 1(a). 

We also detect L1448 C in our map, which is the driving 
source of the high velocity outflow in L1448 (Bachiller et al. 1990). 
The blueshifted CO outflow lobe driven by LI 448 C begins with a 
PA= 159° (Bachiller et al. 1995) (see Figure 1(a)) before being de- 
flected through a total angle of ~ 32° due to a possible collision 
with the core containing L1448 NA and NB (Davis & Smith 1995; 
Curiel et al. 1999). 

HH 7-11. The AMI 1.9 cm data presented here (see Appendix A) 
shows emission extended in a direction almost perpendicular to the 
direction of the HH 7-11 outflow, but consistent with the distribu- 
tion of VLA objects of Rodriguez et al. (1999), almost all of which 
are associated with YSOs. 
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Table 1. Sample selection. Column [1] contains the source name; [2] the protostellar evolutionary class as defined by our criteria described in §2, [3] the Right 
Ascension in units of hours, minutes and seconds; [4] the Declination in units of degrees, minutes, and seconds; [5] the date the observation was made, [6] 
AMI flux calibrator; [7] AMI phase calibrator; [8] rms noise measured from recovered map; [9] major axis of AMI synthesized beam; [10] minor axis of AMI 
synthesized beam. 



Source 


Class 


RA 


Dec 


Date 


1° 


2° 


f'rms 


'^maj 


^min 






(J2000) 


(J2000) 


(dd-mm-yyyy) 






(/ijybeam"') 


(arcsec) 


(arcsec) 


L1448 IRS 3 





03 25 36.49 


+30 45 22.0 


23-09-2010 


3C48 


J0324+3410 


24 


41.8 


27.4 


HH7-11 


I 


03 29 03.73 


+31 16 03.8 


24-09-2010 


3C48 


J0324+3410 


26 


40.2 


28.2 


L1551 IRS 5 


I 


04 31 34.15 


+18 08 04.8 


30-09-2010 


3C48 


J0428+1732 


24 


50.9 


27.2 


L1527 





04 39 53.87 


+26 03 09.9 


25-09-2010 


3C48 


J0440+2728 


19 


44.3 


28.4 


HH1-2MMS1 





05 36 22.85 


-06 46 06.6 


12-10-2010 


3C48 


J0541-0541 


61 


70.5 


24.0 


HH 26 IR 


I 


05 46 03.90 


-00 14 52.5 


27-08-2011 


3C48 


J0558+0044 


34 


78.0 


28.5 


HHlll 





05 51 46.25 


+02 48 29.7 


06-10-2010 


3C48 


J0552+0313 


31 


68.7 


27.4 


NGC 2264 G 





06 41 11.05 


+09 55 59.2 


04-10-2010 


3C48 


J0654+1004 


33 


54.5 


26.9 


Serpens MMS 1 





18 29 49.79 


+01 15 20.8 


11-10-2010 


3C48 


J1824+0119 


36 


66.5 


25.7 


L723 





19 17 53.67 


+19 12 19.6 


30-09-2010 


3C48 


J1914+1636 


23 


48.1 


26.3 


L1251 A 





22 35 24.95 


+75 17 11.4 


25-09-2010 


3C147 


J2236+7322 


25 


37.6 


27.8 



Table 2. Additional detected sources. Column [1] contains the additional source number, [2] the field the source is located in, [3] - [4] the Right Ascension 
and Declination of the peak flux value, [5] the AMI 16 GHz combined channel flux density, [6] the spectral index oIami calculated from the AMI fluxes, [6] 
the associated source, and [7] the protostellar evolutionary class where applicable or object type. 



AMI 


Field 


RA 


Dec 


5l6GHz 


«AMI 


Associated 


Class 


Source 




(J2000) 


(J2000) 


(mJy) 




Source 




1 


L1448 


03 25 38.9 


+30 44 02 


0.76 ±0.09 


-0.67 ±0.98 


L1448 C 





2 


HH7-11 


03 28 55.5 


+31 14 35 


0.39 ±0.10 


1.84±0.16 


NGC 1333 IRAS 2A 





3 


HH7-11 


03 28 57.3 


+31 14 16 


0.50 ±0.10 


1.65 ±0.34 


NGC 1333 IRAS 2B 


I 


4 


L1551 


04 3144.3 


+18 08 29 


0.74 ±0.09 


0.61 ±0.75 


L1551 NE 





5 


L1527 


04 39 49.4 


+26 05 25 


0.96 ±0.09 


-0.10 ±0.92 


HH192VLA2 


ext? 


6 


LI 527 


04 39 59.8 


+26 02 06 


0.88 ±0.09 


-1.25 ±0.58 


HH 192 VLA 3 


ext? 


7 


HH 1-2 


05 36 25.7 


-06 47 16 


1.00 ±0.26 


0.10 ±1.10 


HH2 


HH 


8 


HH 26 IR 


05 46 07.4 


-00 13 42 


1.23±0.10 


1.82±0.13 


HH 25 MMS 





9 


HHlll 


05 5145.5 


+02 50 06 


0.70±0.15 


-0.23 ±1.10 


- 


ext? 


10 


L1251 


22 35 20.5 


+75 19 44 


0.71 ±0.09 


-1.20±0.68 


- 


ext? 



HH 7-11 SVS 13 is a designated Class I object in NGC 1333 
within the Perseus molecular cloud, and is generally believed to be 
the driving source of the chain of HH 7-11 knots to the southeast. 
However, using high-resolution VLA data Rodriguez et al. (1997) 
propose that VLA 3, detected ~ 6" southwest of SVS 13 (VLA 4), 
is instead the driving source of the HH 7-11 outflow. VLA 3 is 
better aligned with the HH 7-11 knots and the elongation of the 
source is approximately along the axis of the HH and molecular 
flows. Rodriguez et al. (1997) suggest that SVS 13 is the driving 
source of a much smaller, nearby chain of HH objects to the east. 

We also detect NGC 1333 IRAS 2A and 2B in our map, which 
were first detected by Jennings et al. (1987) as IRAS 2, and subse- 
quently by Sandell et al. (1994) who detected a quadrupolar out- 
flow centered on the source suggesting that it might be a binary. 
The two sources, IRAS 2A and 2B, were detected at 3.6 and 6 cm 
by Rodriguez et al. (1999) (VLA 7 and 10 respectively) who sug- 
gest that the quadrupolar outflow may originate from IRAS 2A and 
which therefore itself may be a binary. 

L1551. L1551 IRS 5 is a deeply embedded multiple protostellar 
Class I source in Taurus that drives a binary jet (Rodriguez et al. 
2003a). At 1.9 cm (see Appendix A) this source appears point- like 
and it has been suggested that the emission for this source near 



1 cm arises from nearly equal contributions of dust emission from 
the discs and free-free emission from the ionised ejecta (Rodriguez 
1998). 

Rodriguez et al. (2003b) have shown that the emission at cen- 
timetre wavelengths is produced by a pair of closely aligned (within 
~ 12°) bipolar ionised jets, which themselves are aligned with the 
larger scale bipolar molecular outflow. Lim & Takakuwa (2006) de- 
termined that L1551 IRS 5 is a triple protostellar system by identi- 
fying two main (N and S) components, separated by ~ 0.3", each 
comprised of a circumstellar dust disc and bipolar ionised jet, and 
a third component lying ~ 0.09" southeast of the N component as 
another circumstellar dust disc. 

We detect, but do not resolve, the Class source LI 551 NE 
(Moriarty-Schieven et al. 1995) to the east of IRS 5 at 1.9 cm. 
L1551 NE has a molecular outflow and is likely to drive the HH 
objects 28, 29, and 454 (Devine et al. 1999). Radio continuum ob- 
servations at 3.5 cm identified L1551 NE as a binary source with 
separation 0.5" (Reipurth et al. 2002). 

L1527. We detect L1527 (IRAS 04368+2557) at 1.9 cm as a point- 
source (see Appendix A). This object is an embedded Class 
(Melis et al. 2011) protostar in the Taurus molecular cloud with 
a compact accretion disc and a ~ 0.2" binary companion (Loinard 
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Table 3. Source information. Column [1] contains the source name, [2] the molecular cloud association, [3] the distance to the source, [4] distance reference, 
[5] bolometric temperature, [6] bolometric luminosity, [7] reference for Tbol ^ind Lboi, [8] envelope mass calculated using the 850 /im flux from the Fundamental 
Map Object Catalogue (Di Francesco et al. 2008), [9] outflow force of the source, and [10] reference for the outflow force of the source. 



Source 


Cloud 


d 


d 


Tbol 


Lboi 


Ref 


Mcnv 


^out 


.Pout 




Association 


(pc) 


Ref 


(K) 


(Lq) 




(Me) 


(Mg' kms"' yr"') 


Ref. 


L1448 IRS 3 


Perseus 


250 


1 


53 


13.7 


6 


16.22 


49.9x10-2 ±5.29 


12 


HH7-11 


Perseus 


250 


1 


180 


18.0 


6 


10.75 


1.53±18.0 


12 


L1551 IRS 5 


Taurus 


140 


2 


75 


15.3 


7 


6.22 


9.63 


13 


L1527 


Taurus 


140 


2 


59 


2.0 


8 


5.89 


0.79±0.13 


17 


HH1-2MMS 1 


Orion 


420 


3 


- 


23 


9 


42.46 


- 


- 


HH 26 IR 


Orion 


420 


3 


- 


25.1 


10 


4.68 


36.1 


13 


HHlll 


Orion 


420 


3 


49 


19.2 


11 


13.90 


4.57 


13 


NGC 2264 G 


MonOBl 


800 


4 


25 


12 


4 


11.59 


300 


13 


Serpens MMS 1 


Serpens 


260 


14 


56 


17.3" 


15 


28.17 


2.8 


16 


L723 


- 


300 


4 


50 


3 


4 


5.14 


26 


13 


L1251 A 


Cepheus Flare 


320 


5 


- 


2.24 


5 


14.83 


8.6 


13 


L1448 C 


Perseus 


250 


1 


49 


3.9 


6 


5.1 


16.3x10-2 ±1.9 


12 


NGC 1333 IRAS 2A 


Perseus 


250 


1 


58 


13.9 


6 


2.81* 


56.8x10-2 ±5.6 


12 


NGC 1333 IRAS 2B 


Perseus 


250 


1 


106 


5.3 


15 


0.74* 


- 


- 


L1551 NE 


Taurus 


140 


2 


56 


4.2 


11 


5.6 


- 


- 


HH 25 MMS 


Orion 


420 


3 


34 


5.2 


4 


43.96 


- 


- 



References: (1) Dunham et al. 2008; (2) Kenyon et al. 1994; (3) Menten et al. 2007; (4) Andre et al. 2000; (5) Kun et al. 2009; (6) HatcheU 

et al. 2007b; (7) van Kempen et al. 2009; (8) J0rgensen et al. 2009; (9) W. J. Fischer et al. 2010; (10) Nisini et al. 2002; (11) Froebrich 2005; 

(12) Hatchell et al. 2007a; (13) Anglada 1995; (14) Straizys et al. 1996; (15) Enoch et al. 2009b; (16) S. R Graves et al. 2010; (17) Bontemps 

et al. 1996. 

" Enoch et al. 2009a 

* Enoch et al. 2009b 

Table 5. Integrated flux densities. Column [1] contains the source name, [2] the combined-channel 16 GHz flux density for that source, [3]-[7] contain the 
AMI integrated flux densities for each source from AMI channels 4-8, and [8] contains the spectral index cKami calculated from the AMI fluxes. 





5|6GHz 




AMI Channel Number 






Source 


4 


5 


6 


7 


8 


«AMI 






[14.62 GHz] 


[15.37 GHz] 


[16.12GHz] 


[16.87 GHz] 


[17.62 GHz] 






(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 




LI 448 IRS 3 


2.39±0.16 


2.40 ±0.25 


2.34 ±0.26 


2.41 ±0.37 


2.26 ±0.24 


2.62 ±0.37 


0.20 ±0.73 


HH7-11 


3.63 ±0.24 


3.43 ±0.20 


3.10±0.18 


3.91 ±0.22 


3.43±0.19 


3.83 ±0.41 


0.70 ±0.41 


L1551 IRS 5 


4.97 ±0.27 


4.74 ±0.25 


4.37 ±0.24 


4.79 ±0.28 


4.82 ±0.29 


5.48 ±0.57 


0.79 ±0.40 


L1527 


1.33±0.10 


1.04±0.11 


1.04 ±0.14 


1.12±0.12 


1.20±0.16 


1.38±0.19 


1.15 ±0.57 


HH 1-2 MMS 1 


1.70 ±0.27 


1.30±0.28 


1.53 ±0.30 


1.34±0.18 


1.33±0.14 


- 


-0.17±1.02 


HH 26 IR 


0.39 ±0.08 


- 


- 


- 


- 


- 


- 


HHlll 


2.49 ±0.20 


3.11±0.53 


2.68±0.28 


2.39±0.26 


2.86±0.30 


2.51 ±0.40 


-0.40 ±0.85 


NGC 2264 G 


0.30 ±0.22 


- 


- 


- 


- 


- 


- 


Serpens MMS 1 


7.32±0.41 


8.21 ±0.77 


7.66±0.55 


7.04 ±0.43 


6.77 ±0.44 


7.30 ±0.83 


-0.73 ±0.56 


L723 


0.58 ±0.09 


0.58 ±0.09 


0.52 ±0.06 


0.65 ±0.08 


0.66 ±0.07 


0.52 ±0.11 


0.42 ±0.82 


L1251 A 


1.05 ±0.12 


0.86 ±0.18 


0.98 ±0.14 


1.00 ±0.20 


0.76±0.16 


1.10 ±0.25 


0.06 ±0.99 



et al. 2002). The HH 192 outflow associated with this object has an 
east- west oiientation (Gomez et al. 1997) with an edge-on disc per- 
pendicular to the outflow axis (Tobin et al. 2010) that is infalling 
onto the central source (Ohashi et al. 1997). AIVII Consortium: 
Scaife et al. (2012a) found an opacity index jS = 0.32 ± 0.04 for 
this source, consistent with the presence of large dust grains, or 
pebbles, in the disc. 

We also detect HH 192 VLA 2 and 3 (Rodriguez & Reipurth 
1998) at 1.9 cm; there are no other references to these objects in 
the literature. We find a spectral index of 0!ami = —0. 10 for VLA 2 
consistent with optically thin free-free emission and aAMi = — 1 .25 



for VLA 3 which is inconsistent with free-free emission and sug- 
gests that it could be extragalactic. For the purposes of this work, 
VLA 2 and 3 are considered to be extragalactic sources (see §6) 
although the emission from VLA 2 could be consistent with an HII 
region. 

HH 1-2. At 1.9 cm, we detect the HH 1-2 outflow driving source 
as a point source (see Appendix A). HH 1 and 2 are considered 
to be the prototypes of the HH phenomena, displaying two bright 
bow shocks which expand as a part of a bipolar outflow moving 
away from the central source that drives a finely collimated jet 
(Rodriguez et al. 2000). The Class driving source of HH 1 and 



6 Ainsworth et al. 







-' ' ' 


1 1 1 1 1 




- 




+30°48' 




\^ 




- 


o 
o 


+46' 


/ 
/ 


\ 


\ 
\ 


; 


o 




1 
1 

; \ 


M ' 


/ 


- 


01 
Q 


+44' 


- \ / ( 

: \ / 






- 




+42' 


^ 


,,,,!,, 







Right Ascension (J2000) 

(a) 




Figure 1. L1448. (a) The AMI 16 GHz combined channel map uncorrected 
for the primary beam response for L1448. The positions of known Class 
and Class I objects are indicated as crosses (+) and stars (*) respectively, 
and are from Hatchell et al. (2007b). The solid black lines indicate out- 
flow axes for L1448 NB (northern axis) and L1448 C (southern axis) from 
Wolf-Chase et al. (2000), details in §5. We plot the half power point of the 
primary beam as a solid circle (ai 6arcmin at 16 GHz) and the FWHM of 
the PSF as the filled ellipse in the bottom left corner (see Table 1). Con- 
tours at 5, 10, 15, 20c7niis etc, where CTms = 24^Jybeam"'. (b) The SED 
for L1448 IRS 3 with the observed radio spectra over the AMI frequency 
channels 3-8 combined with flux densities from the literature. 



2, VLA 1, was detected as a centimetre continuum radio source by 
Pravdo et al. (1985), and subsequent observations have shown that 
the source powers a small thermal radio jet along the axis of the HH 
flow (Rodriguez et al. 1990). Another source, VLA 2, was found 3" 
from VLA 1, producing the HH 144 flow which is at a large angle 
to the HH 1 flow axis. The AMI does not resolve the two sources 
VLA 1 and 2 and the combined thermal emission is approximately 
point-like. 

The Herbig-Haro objects HH 1 and HH 2 themselves have 
been previously detected at centimetre wavelengths (e.g. Rodriguez 
et al. 1990, 2000) and we detect HH 2 but we do not detect HH 1 
at 1.9cm. Rodriguez et al. (1990) find a flux density of 0.78 mjy at 
2 cm for HH 1, however we constrain an upper flux density limit 
of 0.34 mJy, suggesting that this source is variable. Radio contin- 
uum emission from other HH objects has been found to be time 
variable (Rodriguez et al. 2008a; Carrasco-Gonzalez et al. 2010b). 



Rodriguez et al. (1990) find flux densities of 1.75 mJy for the cen- 
tral source and 0.98 mJy for HH 2, for which our data is consistent. 
HH 26 IR. HH 26 IR is a weak source at 1.9 cm (see Appendix A), 
a Class I protostar (Davis et al. 1997) located within L1630 in Orion 
that drives a molecular outflow. 

LI 630 is rich in star formation activity and has a complex 
morphology. Two embedded Class protostars, HH 24MMS and 
HH 25MMS, are also in this star formation region and both have 
compact jets, while the outflow of HH 26 IR drives a more ex- 
tended outflow (Gibb & Heaton 1993). H2 knots HH 26A, B, C, 
and D exist along the axis of the HH 26 IR outflow which is ap- 
proximately orthogonal to the HH 25 outflow (Davis et al. 1997). 
We also detect HH 25 MMS within the AMI primary beam to the 
northeast of HH 26 IR. 

HH 111. At 1.9cm, HH 111 (see Appendix A) is significantly 
extended along an approximate north-south direction, perpen- 
dicular to the known outflow direction. The driving source is 
IRAS 05491-1-0247, a suspected Class I binary source (Yang et al. 
1997) that drives a bipolar molecular outflow. The outflow was first 
discovered in the optical by Reipurth (1989). However, based on 
the classification criterion defined in §2, this source has a bolomet- 
ric temperature consistent with a Class object and this is discussed 
further in §6. The jet is well-collimated, aligned approximately in 
the E-W direction, and has knots with velocities ranging from 300 
to 600 km s^ ' for different parts of the jet complex (Reipurth et al. 
1992). It was found that the optical jet is part of a giant HH complex 
extending over 7.7 pc (Reipurth et al. 1997). 

Near-infrared observations (Gredel & Reipurth 1994) revealed 
a second bipolar flow, HH 121, emerging from around the same 
position as the optical outflow and is aligned approximately in the 
N-S direction, consistent with the extension observed here. Such 
a quadrupolar outflow is indicative of a binary companion. 3.6 cm 
(Reipurth et al. 1999) and 7 mm (Rodriguez et al. 2008b) VLA ob- 
servations suggest a common origin within ~ 0. 1". It has been sug- 
gested that the most viable interpretations for the structure are that 
the observations are of two orthogonal discs around separate proto- 
stars or a disc with a perpendicular jet (Rodriguez et al. 2008b). 
NGC 2264. We detect emission from NGC 2264 G at 9c7rms in the 
AMI combined-channel map (see Appendix A). The NGC 2264 G 
molecular outflow was one of nine discovered in a CO survey of 
the Monoceros OBI cloud conducted by Margulis & Lada (1986) 
and was found to have very high velocity emission and a well- 
collimated, bipolar morphology (Margulis et al. 1988). It is cen- 
tered on IRAS 06384-1-0958 and extends ~1.6pc in an approxi- 
mately E-W direction (Caratti o Garatti et al. 2006). Lada & Fich 
(1996) show that the outflow is S-shaped and symmetric with re- 
spect to the driving source, VLA 2, which was detected by Gomez 
et al. (1994) using VLA ammonia observations. They obtained a 
spectral index for the outflow of a ~ 0.3 ± 0.2 which is con- 
sistent with optically thin free-free emission of a thermal jet, and 
their 3.6 cm VLA images show VLA 2 to be elongated in the direc- 
tion of the outflow. Subsequent observations of VLA 2 by Ward- 
Thompson et al. (1995) using the James Clerk Maxwell Telescope 
(JCMT) obtained a ratio of submillimeter and bolometric luminos- 
ity indicative of a Class protostar. 

Serpens. At 1.9 cm, AMI does not resolve the triple radio source 
SMM 1 but measures the integrated radio emission. The emission 
is oriented in an approximate N-S direction (see Appendix A), de- 
spite the outflow orientation of NW-SE which could be due to the 
presence of multiple objects in this region. The Serpens molecular 
cloud is a nearby star formation region containing clusters of low- 
mass protostars as well as several HH objects and outflows. The 
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most luminous and most deeply embedded object in the Serpens 
cloud core is the triple radio source SMM 1, comprised of a cen- 
tral source, a NW component, and a SE component, and is a known 
Class source associated with a highly coUimated radio continuum 
jet (Rodriguez et al. 1989b). 



Curiel et al. (1993) presented high-sensitivity, multifrequency 
VLA radio continuum observations of the triple radio source, show- 
ing a one-sided radio jet morphology, but did not detect a counter 
jet. These authors show that the central source is extended at 
3.6 cm, but the difference between the direction of the major axis 
and the orientation of the radio jet appear to be significant, suggest- 
ing that the central source is precessing or nutating. They find that 
the central source has a spectral index a ~ 0.15 ±0.09, consis- 
tent with the value ofa~0.1±0.1 obtained by Rodriguez et al. 
(1989b), while the NW and SE components have spectral indicies 
a ~ -0.05 ±0.05 and a ~ -0.30±0.04 respectively. These results 
suggest a thermal origin for the emission associated with the cen- 
tral source and the NW component, and a borderline thermal/non- 
thermal origin for the SE component. AMI Consortium: Scaife 
et al. (2012b) present evidence for radio variability in SMM 1, as 
previously suggested by Enoch et al. (2009a) who detect a high 
mass disc around the source. Here we adopt a bolometric luminos- 
ity of Lboi = 17.3 Lq (Enoch et al. 2009a), however we note that 
significantly higher luminosities have been reported for this object 
(e.g. 46 L© from Hurt & Barsony 1996). 

L723. L723, also known as IRAS 19156-1-1906, appears point-like 
at 1.9 cm (see Appendix A). L723 is a designated low-mass. Class 
object (Lee et al. 2002). The associated outflow consists of a pair 
of bipolar lobes aligned along the east-west direction and another 
pair of bipolar lobes aligned roughly in the north-south direction, 
with IRAS 19156-1-1906 located at their common centre (Lee et al. 
2002). VLA observations at 3.6 cm reveal two sources, VLA 1 and 
VLA 2, toward the centre of the outflow (Anglada et al. 1991), al- 
though VLA 1 is likely a background source (Anglada et al. 1996; 
Girart et al. 1997). VLA 2 has a jet-like morphology and partially 
optically thick free-free emission, characteristic of a thermal radio 
jet, and was first identified as the powering source of the EW pair 
of molecular lobes by Anglada et al. (1996). Recent VLA obser- 
vations by Carrasco-Gonzalez et al. (2008) at 3.6 cm and 7 mm re- 
solve VLA 2 into several components and the two brightest sources 
at 3.6cm, VLA 2A and 2B, are separated by Ri 0.29". 



L1251. The emission at 1.9 cm of L1251 A is detected as a point- 
source (see Appendix A). LI 251 A is a star- forming region toward 
Cepheus, and is centered on the bright source IRAS 22343-1-7501 
which is associated with a large-scale (~ 10') molecular outflow 
(Sato & Fukui 1989; Schwartz et al. 1988). VLA observations 
(Meehan et al. 1998; Beltran et al. 2001; Reipurth et al. 2004a) 
show several radio sources (VLA 6, 7, and 10) within this region 
that could possibly drive the molecular outflow. Reipurth et al. 
(2004a) show two sources with clear extension, VLA 7 and 10, 
and suggest that VLA 10 could drive the molecular outflow found 
by Sato & Fukui (1989), which extends for 10' in the approximate 
NE-SW direction. However, VLA 7 is consistent with the outflow 
found by Schwartz et al. (1988) that has an outflow direction in an 
almost perpendicular direction, and the emission we detect is con- 
sistent with this. We do not resolve these sources, but measure the 
integrated radio emission. 



Table 6. Summary of detection statistics. The original sample includes only 
objects listed in Table 1. The extended sample includes all objects listed in 
Tables 1 and 2, inclusive of the original sample. 



Class 


Present 


Detected 


per cent 


Original Sample: 
Class 
Class I 


8 
3 


8 
3 


100 
100 


Extended Sample: 
Class 
Class I 


15 

7 


12 

4 


80 

57 



6 DISCUSSION 

6.1 Detections and non-detections 

We detect 100 per cent of the objects in the target sample with AMI 
at 16 GHz, which includes eight Class and three Class I protostars 
based on the classification scheme described in §2. We also detect 
five additional Class and one additional Class I sources within the 
AMI primary beam of these fields, all of which are Class sources 
driving outflows except for NGC 1333 IRAS 2B which is Class I. 

In the LI 448 field there is one Class source (31; Hate hell 
et al. 2007b) that we do not detect, in the HH 7-11 field there is one 
Class (62; Hatchell et al. 2007b) and one Class I (50; Hatchell 
et al. 2007b) that we do not detect, and in the Serpens field there is 
one Class (VLA 5; AMI Consortium: Scaife et al. 2012b) and two 
Class I (VLA 9, 16; AMI Consortium: Scaife et al. 2012b) sources 
that we do not detect. 

We identify AMI 7 as the radio counterpart to HH 2, but we 
do not detect HH I at 1 .9 cm even though it has been previously de- 
tected at centimetre wavelengths. We also do not detect any HH 192 
objects in the LI 527 field which lie within the AMI primary beam. 

Detection statistics for the original and extended samples di- 
vided by evolutionary class are listed in Table 6. It is not surprising 
that we should detect 100 per cent of the target sample as all sources 
have been previously detected at centimetre wavelengths and have 
known radio jets. Detection rates for Class in the extended sam- 
ple is typical compared to AMI Consortium: Scaife et al. (2011b, 
2012b) but high for Class I. This is likely due to the small sample 
size or the fact that a majority of the extended sample consists of 
the original sample which was chosen on the basis that these ob- 
jects have been detected at these wavelengths before. The original 
sample is biased, and for these reasons it is difficult to draw strong 
statistical conclusions from this sample. 



6.2 Expected contamination by extragalactic radio sources 

At 16 GHz we expect a certain number of extragalactic radio 
sources to be seen within each of our fields. Following AMI Con- 
sortium: Scaife et al. (201Ia,b, 2012b), to quantify this num- 
ber we use the 15 GHz source counts model from de Zotti et al. 
(2005) scaled to the IOC survey source counts (AMI Consortium: 
Davies et al. 2011). The average rms noise from our datasets is 
^ 30/xJy beam^' and from this model we predict that we should 
see 0.043 sources arcmin^ , or b^ 2 radio sources within a 6 arcmin 
FWHM primary beam above a 5c7nns flux density of 150 /iJy. Ac- 
counting for the radial attenuation of the primary beam, we would 
therefore expect to see ~ 10± 3 extragalactic radio sources within 
our target fields. Making the assumption that all sources which can- 
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not be identified witli a previously laiown protostellar object are 
extragalactic we find 4 radio sources, lower than predicted but con- 
sistent at the 2(T level. The low number of extragalactic sources 
detected in these fields compared to the model may represent an 
over-estimation from the 15 GHz source counts, which are extrap- 
olated from a completeness limit of 0.5 mjy and is high compared 
with our detection limit of 0. 15mJy. One break in the observed 
source counts is already known at mJy flux density levels (AMI 
Consortium: Davies et al. 2011) and it is possible that a further 
break exists below 0.5 mJy. 



6.3 Spectral Energy Distributions 

An extensive literature search was conducted for unresolved, inte- 
grated flux densities to include in the spectral energy distributions. 
It should be noted that Wendker (1995) was a useful reference. 
High resolution data that were highly discrepant due to flux loss or 
data with high uncertainties (e.g. 450 /im data from di Francesco 
et al. 2008) were not included. Spectral energy distributions are 
shown in Figures 1(b), 2, and Appendix B, with maximum likeli- 
hood models (see Table 8) overlaid. The list of archival data used 
in the spectral energy distributions can be found in Appendix C. 
Where uncertainties were not provided, an error of 10 per cent was 
used in the model fittings and this is indicated by a ^ in Appendix C. 
Only data v < 3THz (A > 100/lm) were included in the fit, but 
IRAS data v > 3 THz are included in the plots for illustration. 



6.3.1 Maximum likelihood SED fitting 

Power-law spectral indices, ojami^ were fitted to the AMI chan- 
nel data alone for each object, see Tables 2 and 5, using the 
Markov Chain Monte Carlo based Maximum Likelihood algorithm 
METRO (Hobson & Baldwin 2004). The fit is of the form: 



Sv o= V 



O^AMI 



(1) 



where Sy is the flux density measured at frequency v. 

METRO was then used to fit a combined radio power law, 
with spectral index a', and blackbody model to the larger dataset 
for each source. This fit utilized data at wavelengths longer than 
100 /im and had the form: 



'5'total —S\+S2— Ki 






-K2 



v^ByjTi) 



(2) 



where /3 is the dust opacity index, Sy is the Planck function for 
a dust temperature 7^, Ki is the normalized flux density at Vj = 
16 GHz and K2 is the normalized flux density at V2 = 300 GHz. Fit- 
ting was performed in two scenarios to obtain the maximum like- 
lihood values for the spectral index a' of each source. These two 
scenarios are: 

(i) fixing the dust temperature based on evolutionary class, 
(ii) allowing the dust temperature to vary. 

It can be seen that when v = Vi , Si equals Ki , the normalized flux 
density at 16 GHz and when v = V2, 52 equals K2, the normalized 
flux density at 300 GHz, and we define these parameters as 5""™ 
and S^o™ respectively. We use uniform and separable priors for all 
parameters, with ranges 



6.3.2 The grey body contribution 

When V = Vi = 16 GHz, the free-free component should dominate 
the total flux density, however at 16 GHz there is also expected to 
be a small contribution to the radio flux density of protostars due 
to the long wavelength tail of the thermal dust emission from the 
envelopes around these objects (e.g. AMI Consortium: Scaife et al. 
2012a). 

It is important to separate the greybody component from the 
free-free emission and we calculate this predicted contribution to 
the 16 GHz flux density in the following way: 



■5gb,16 =f{K2,l^,Td) = ^2—p , 

^2^300(^(1) 



with an associated error 
dK2 



3f. 



dTA 



f^Ji 



(4) 



(5) 



n = n„(-2,2)n^(o,3)n7i(5,45) 



(3) 



where v = 16 GHz, V2 = 300 GHz, and dK^, <yp and aj;, are the 
uncertainties on K2, /3 and T^ respectively and listed in Tables 7 
and 8. In all cases, we find that the error, CT5 jg, is dominated by 
the uncertainty on j8. For example, in the case of LI 448: csf ig = 

4.392 ±75.82 ±14.3 Vjy- 

We subtract the value of 5gb.i6 from the measured 16 GHz flux 
density to remove the contribution of thermal dust emission and 
obtain values of only the radio emission due to the outflow com- 
ponent at 16 GHz, 5rad.i6- This is to ensure that the values used do 
not include contributions from the thermal dust tail which varies 
greatly between sources and which might therefore influence any 
conclusions being drawn from the correlations examined in a later 
section of this paper. We compute the radio luminosity S^^^i^d 
using the greybody-subtracted radio flux densities and the distance 
values listed in Table 3 to use in the correlations in §6.6. If the grey- 
body subtracted flux density is not distinct from zero by > 3(7, we 
consider the radio luminosity consistent with zero (see Tables 7 and 
8). 



6.3.3 Comparison of the two scenarios 

In Scenario (i), the algorithm was run by fixing T^ to 15 K for 
Class objects and to 20 K for Class I objects, and allowing the 
other parameters to vary to determine their maximum likelihood 
values. Overall, fixing the dust temperature provided reasonable fits 
to the spectral energy distributions, with the exception of L1527, 
NGC 2264 G, Serpens MMS 1, and L1551 NE where it provided a 
poor fit to the greybody peaks. The results for the maximum like- 
lihood values of a', j8, 5™™ and S^g™ in each case are listed 
in Table 7, along with the dust temperature used [Columns 2-6], 
the greybody contribution at 16 GHz [Column 7], and the radio 
luminosity at 16 GHz with the greybody contribution subtracted 
out [Column 8]. We find a weighted average spectral index of 
a' — 0.20 ±0.41 for the original sample, consistent with free-free 
emission, and a dust opacity index of jS = 1 .75 ± 0.36. Quoted er- 
rors represent the standard deviation of the distribution of values, 
rather than the standard error on the mean. 

In Scenario (ii), we do not fix T^ and allow the algorithm to 
determine the maximum likelihood value for this parameter along 
with a', j3, SJ™™ and S^g™. The results for these parameters along 
with 5gb,i6 and Sj-ad.ie"^^ are listed in Table 8. We again find a 
weighted average spectral index of a' = 0.20 ±0.46 for the origi- 
nal sample consistent with free-free emission. We find an opacity 
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Table 7. Model results for fixed Tj. Column [1] contains the source name, [2] the spectral index, [3] the opacity index, [4] the dust temperature, where 15 K 
was used for Class sources and 20 K for Class I, [5] the normalized flux density at 16 GHz, [6] the normalized flux density at 300 GHz, [7] the predicted 
greybody contribution at 16 GHz, and [8] the radio luminosity measured at 16 GHz with the predicted thermal dust contribution subtracted. 'UC'is used to 
indicate that the parameter is unconstrained. 



Source 


a' 


P 


Ta 


onorm 


cTiorm 
"^300 


Sgb.16 


^rad.16^ 








(K) 


(mjy) 


(Jy) 


(MJy) 


(mJykpc^) 


LI 448 IRS 3 


0.37 ±0.09 


1.48 ±0.05 


15 


1.88 ±0.09 


4.96±0.15 


297.47 ±52.50 


0.131 


HH7-11 


1.12±0.04 


1.92 ±0.02 


20 


3.03 ±0.08 


5.30±0.19 


77.96 ±4.83 


0.222 


L1551IRS5 


0.10 ±0.04 


1.52 ±0.14 


20 


3.90±0.17 


4.47 ±0.43 


214.58± 130.38 


0.093 


L1527 


0.17 ±0.08 


0.96 ±0.10 


15 


0.85 ±0.05 


0.49 ±0.02 


134.94±52.75 


0.023 


HH1-2MMS 1 


0.23 ±0.04 


2.38 ±0.08 


15 


1.49 ±0.07 


0.99 ±0.02 


4.33±1.25 


0.299 


HH 26 IR 


UC 


0.62 ±0.07 


20 


0.04 ±0.06 


0.56 ±0.03 


372.00 ±101.97 


- 


HHlll 


0.88 ±0.07 


1.66 ±0.07 


15 


2. 16 ±0.09 


0.90 ±0.02 


31.95±8.81 


0.434 


NGC 2264 G 


-0.29 ±0.24 


1.44 ±0.39 


15 


0.50±0.10 


0.54 ±0.65 


37.06 ±192.75 


- 


Serpens MMS 1 


-0.06 ±0.02 


1.47 ±0.05 


15 


6.89±0.14 


5.22 ±0.20 


327.71 ±56.62 


0.473 


L723 


-0.26 ±0.09 


1.83 ±0.07 


15 


0.58 ±0.03 


0.93 ±0.05 


20.04±5.15 


0.050 


L1251 A 


1.79±0.15 


2.92 ±0.27 


15 


0.91 ±0.07 


0.32 ±0.07 


0.29±0.51 


0.107 


L1448 C 


1.12±0.12 


2. 10 ±0.04 


15 


0.53 ±0.02 


1.64 ±0.04 


16.25 ±2.41 


0.046 


NGC 1333 IRAS 2A 


2.42 ±0.07 


2.01 ±0.05 


15 


0.43 ±0.02 


1.52±0.10 


19.33 ±3.67 


0.023 


NGC 1333 IRAS 2B 


1.31 ±0.05 


1.95 ±0.09 


20 


1.02 ±0.04 


0.53 ±0.04 


7.16±2.69 


0.031 


L1551 NE 


-0.13±0.31 


1.41 ±0.03 


15 


0.60 ±0.09 


1.83 ±0.02 


136.09 ±11.08 


0.012 


HH 25 MMS 


2.13 ±0.04 


1.68±0.16 


15 


0.97 ±0.07 


0.83±0.10 


27.82 ±22.30 


0.212 



Table 8. Model results for variable Tj. Column [1] contains the source name, [2] the spectral index, [3] the opacity index, [4] the dust temperature, [5] the 
normalized flux density at 16 GHz, [6] the normalized flux density at 300 GHz, [7] the predicted greybody contribution at 16 GHz, and [8] the radio luminosity 
measured at 16 GHz with the predicted thermal dust contribution subtracted. 'UC'is used to indicate that the parameter is unconstrained. 



Source 



(K) 



(mJy) 



(Jy) 



Sgb.16 

(MJy) 



Smd.ied 

(mJy kpc^) 



L1448 IRS 3 


0.46 ±0.09 


1.81±0.12 


10.77±0.95 


2.02 ±0.09 


5.13±0.12 


145.74±77.31 


0.140 


HH7-11 


1.28 ±0.04 


3.00 ±0.08 


12.53 ±0.90 


3.24 ±0.07 


6.21 ±0.27 


4.92 ±1.59 


0.227 


L1551IRS5 


0.09 ±0.03 


1.46 ±0.05 


26.74 ±2.07 


3.87 ±0.0.09 


4.54±0.17 


233.58±43.54 


0.093 


LI 527 


0.20 ±0.08 


0.81±0.11 


44.86 ±1.57 


0.86 ±0.06 


0.44 ±0.03 


138.81 ±60.00 


0.023 


HH 1-2 MMS 1 


0.21 ±0.04 


2.18±0.19 


16.93 ±1.80 


1.46 ±0.07 


1.01 ±0.02 


7.35 ±7.38 


0.299 


HH 26 IR 


UC 


0.66 ±0.16 


17.36±5.25 


0.04 ±0.09 


0.59 ±0.04 


368.59 ±300.38 


- 


HHlll 


0.84 ±0.07 


1.55±0.17 


17.33 ±3.68 


2.11±0.10 


0.90 ±0.02 


41.31 ±35.89 


0.366 


NGC 2264 G 


-0.31±0.11 


1.71±0.13 


18.30±1.88 


0.49 ±0.06 


0.36 ±0.02 


10.02 ±5.59 


0.186 


Serpens 


-0.07 ±0.02 


1.25 ±0.04 


28.54 ±1.54 


6.74 ±0.14 


4.85±0.18 


451.87±71.01 


0.464 


L723 


-0.29 ±0.09 


1.59±0.13 


17.56±1.31 


0.56 ±0.04 


0.92 ±0.04 


37.73 ±20.70 


0.049 


L1251 A 


1.70 ±0.14 


2.27 ±0.41 


20.85 ±4.36 


0.86 ±0.07 


0.38 ±0.06 


1.98±7.84 


0.107 


L1448 C 


1.17±0.15 


2.22 ±0.11 


12.70±1.25 


0.54 ±0.03 


1.69 ±0.04 


12.88 ±5.99 


0.047 


NGC 1333 IRAS 2A 


2.56 ±0.03 


2.76±0.19 


9.85 ±1.10 


0.46 ±0.02 


1.17 ±0.07 


2.22 ±2.20 


0.024 


NGC 1333 IRAS 2B 


1.40 ±0.05 


2.89 ±0.26 


9.55 ±1.90 


1.06 ±0.03 


0.57 ± 0.06 


0.76 ±1.28 


0.031 


L1551NE 


UC 


0.78 ±0.05 


40.64 ±3.29 


0.17 ±0.07 


1.69 ±0.02 


585.75 ±0.05 


- 


HH 25 MMS 


UC 


0.55 ±0.09 


41.27 ±3.65 


0.004 ±0.18 


1.51 ±0.09 


1021.76±352.26 


- 



index of /3 = 1.52±0.56 and a dust temperature T^ = 18.46±9.77 
for the original sample. 

In the first scenario, with fixed T^, 80 per cent of the sources 
in the original sample for which the spectral indices could be con- 
strained had maximum likelihood values of 2 > a' > — 0. 1 , consis- 
tent with free-free emission. In the second scenario, with Tj vary- 
ing, again 80 per cent of the sources in the original sample had 
a' consistent with free-free emission. In both cases, NGC 2264 
and L723 have spectral indices that are inconsistent with free-free 
emission and a' < —0.1, suggesting non-thermal emission. The 
one source for which an accurate value of a' is unconstrained in 



the original sample is HH 26 IR, due to insufficient low-frequency 
data. Unconstrained parameters are indicated in Tables 7 and 8 by 
'UC. Since a' is unconstrained for HH 26 IR, the algorithm at- 
tributes all the flux density at 16 GHz to thermal dust emission, 
resulting in a value of S^-^i\i(,d ~ when subtracted from the AMI 
flux density. Consequently it does not appear in the correlations in 
§6.6. 

A spectral index of a > 2 is generally considered un-physical 
for free-free emission, and is attributed to dust emission. In the 
cases where the likelihood distribution pushed the upper limit of 
the prior range, we increased the prior range to Ha (—2,3). In the 
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extended sample, NGC 1333 IRAS 2A has a' > 2 in both model 
scenarios, suggesting that the lower frequency data is an extension 
of the greybody emission and that there is no free-free component. 
This will not be reflected accurately in the calculation of the grey- 
body emission at 16 GHz since the algorithm constrains a separate 
a' component, and this will in turn affect the radio luminosity cal- 
culated for this source and the correlations discussed in §6.6. When 
Tj is fixed (Scenario (i)), HH 25 MMS also has a' > 2, suggesting 
that it is purely greybody emission. This will have the same conse- 
quences as for NGC 1333 IRAS 2A. However, a' is unconstrained 
for HH 25 MMS when T^ is allowed to vary (Scenario (ii)) and this 
will have the reverse affect on the radio luminosity. Examining the 
SED plot (see Appendix B), it is possible that HH 25 MMS has a 
free-free component based on one or two data points, but the model 
is unable to constrain it without more low-frequency data. 

For LI551 NE, a' is constrained in Scenario (i) but not in 
Scenario (ii). Consequently, the calculated radio luminosities dif- 
fer significantly between the two scenarios and this will affect the 
correlations in §6.6. Specifically, Sj-ad,l6<^^ 's a factor of 4.3 higher 
when Td is allowed to vary than when it is fixed to 15 K (see Ta- 
bles 7 and 8). When Tj is fixed for this source, the peak of the 
greybody is constrained by the given dust temperature (and does 
not fit the higher frequency data), allowing the algorithm to deter- 
mine a' using the two low-frequency data points (see Figure 2(a)). 
However, when T^ is allowed to vary, the AMI data point at 16 GHz 
is fit as part of the greybody emission and a' is left unconstrained, 
as there is only one remaining low-frequency data point (see Fig- 
ure 2(b)). This results in a more significant greybody contribution at 
16 GHz. Therefore the radio luminosity will be significantly lower 
when iSgb.ie is subtracted in this scenario compared to that when the 
dust temperature is fixed, and this will drastically change its posi- 
tion on the correlation plots between the two scenarios. The same 
situation is also found in the case of HH 25 MMS. 





Figure 2. The spectral energy distribution fits for LI 551 NE when (a) 
Ti = 15 K, and (b) when Tj is allowed to vary. In the first case, the method 
determines a' for the two low-frequency data points, but in the second case 
the AMI data point is fit as greybody emission and a' is unconstrained. 



6.4 Dust Temperatures 

In Scenario (ii), where T^ varies, our results for Tj range between 
9^5 K. Values in the literature for the dust temperatures associ- 
ated with Class and I sources are found to range between 10- 
40 K (Reipurth et al. 1993; Dent et al. 1998; Shirley et al. 2002; 
Young et al. 2003; Hatchell et al. 2007b). Dent et al. (1998) find 
that their sample of Class protostars can be fitted with optically 
thin emission, a j3 index between 1-1.5, and a dust temperature be- 
tween 20-30 K. These authors derive for their Class I sources /3 
between 1-2 and a dust temperature of 30 K, of which our results 
are consistent. 

For Class sources, we find weighted average values of j3 = 
1.35 ± 0.29 and 7^" = 20. 1 ± 10.89, and for Class I sources we find 
jS = 1.81 ±0.74 and 7d = 14.85 ±5. 15. These results are reversed 
from what we expect based on the first SED scenario, but this is 
likely to be due to the low number of Class I sources in our sample 
compared to Class 0. AMI Consortium: Scaife et al. (2012a) find 
jS = 0.3 for LI 527 by fitting data from 5 GHz to 350 GHz jointly 
using the sum of two power-laws, which is slightly shallower than 
the value of jS = 0.81 ± 0. 1 1 found in this work. However, both are 
consistent with the value of j8 ^ 1 , indicating the presence of large 
dust grains. 



6.5 Radio Spectral Indices 

Typically, the radio emission detected at centimetre wavelengths 
for Class and I protostars is observed to possess a flat or rising 



spectrum, indicating that it occurs as a consequence of free-free ra- 
diation from ionised gas. Unresolved sources can exhibit partially 
opaque spectra with spectral indices in the range —0.1 ^ a ^ 2, 
where a value of —0.1 indicates optically thin free-free emission 
and a value of 2 indicates optically thick. For the case of a stan- 
dard, canonical jet with constant opening angle, the spectral index 
a = 0.6. Shallower spectral indices (a < 0.6) are expected for unre- 
solved, well-collimated, ionised outflows, and steeper spectra still 
(a > 0.9) are expected for fully ionised, adiabatic jets (Reynolds 
1986). In this section we compare our results for aAMi (Tables 2 
and 5) and a' from the spectral energy distribution scenario where 
Tj is allowed to vary (Table 8). 

Across the AMI band, 7/9 (78 per cent) of the target sources, 
where the spectral index (Zami could be constrained, are within the 
range —0. 1 $: 0!ami ^ 2 indicating free-free emission (see Table 5). 
When including the additional protostellar sources detected in these 
fields, 1 1/14 (79 per cent) of the sources have oJami within this 
range. The three sources that lie outside this range, HH 111, Ser- 
pens and L1448 C, have ttAMl < —0.1, typically indicating that 
another emission mechanism is at work. However, spectral indices 
inconsistent with free-free emission have been found previously 
for Serpens MMS 1 (Rodriguez et al. 1989b; Curiel et al. 1993). 
The weighted average spectral index of the original sample calcu- 
lated from the AMI fluxes is OamT = 0.42 ± 0.59, consistent with 
a canonical jet. The flux densities for HH 26 IR and NGC 2264 
could only be extracted from the 16 GHz combined-channel map. 
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so it was not possible to measure CHami for tiiese sources. For the 
extended sample, Oami = 1 .52 ± 0.64, still consistent with free-free 
emission. The large uncertainties of oJami in Table 5 are due to the 
short frequency coverage of AMT 

For comparison, 8/10 (80 per cent) of the original sam- 
ple have — 0. 1 ^ o:' ^ 2 when including archival data (again, 
excluding HH 26 IR, see Table 8). For the extended sample, 
10/13 (77 per cent) of the sources (excluding L1551 NE and 
HH 25 MMS) have a' within this range. The values for a' that 
lie outside this range are NGC 2264 and L723 (a' < -0.1) and 
NGC 1333 IRAS 2A (a' > 2). The original sample has a weighted 
average spectral index a' = 0.20 ±0.46 and the extended sample 
has a' = 0.37 ±0.63, both consistent with free-free emission and a 
well-collimated flow. 

In the case of NGC 2264 G, there are two known radio con- 
tinuum sources in this outflow region: VLA 1, which was proposed 
as the driving source by Rodriguez & Curiel (1989), and VLA 2, 
which is suggested to be the driving source source by Gomez et al. 
(1994). Gomez et al. (1994) find a spectral index a ~ -1.3 ±0.4 
for VLA 1 and a significantly lower flux density than Rodriguez 
& Curiel (1989), suggesting that it appears to be highly variable 
and an extragalactic nature cannot be ruled out. For VLA 2, these 
authors find a spectral index a ~ 0.3 ± 0.2, consistent with free- 
free emission. The negative spectral index a' = —0.31 ±0.11 ob- 
tained here for NGC 2264 G is likely contaminated by VLA 1 as 
the sources are unresolved by AMI. Additional high resolution ob- 
servations are needed to confirm the non-thermal component. 

L723 is in a similar situation to NGC 2264 G, where there is 
a nearby source, VLA 1, with spectral index a ^ —0.1 that could 
be contaminating the results in the SED for the prospective driving 
source, VLA 2. Anglada et al. (1996) find a spectral index a = 
—0.3 ± 0.2 for VLA 1, and discuss the possibility that the source is 
a radio-emitting, optically obscured T Tauri star. These authors find 
a spectral index a = 0.5 ± 0.2 for VLA 2, consistent with free-free 
emission, and the results of Carrasco-Gonzalez et al. (2008) for the 
spectral index of VLA 2 are also indicative of free-free emission. 
We find Kami = 0.42 ± 0.82 for L723, consistent with free-free 
emission from a canonical jet even though we do not resolve these 
sources. The negative spectral index a' = —0.29 ±0.09 we obtain 
is likely affected by both VLA 1 and insufficient low-frequency 
data to better constrain a'. 

From the data presented here, LI 527 is found to have a 
spectral index chami = 1-15 ±0.57, which agrees with the value 
of chami = 1.17 ±0.42 found by AMI Consortium: Scaife et al. 
(2012a). A spectral index of a > 0.8 can be easily produced by 
velocity gradients or the recombination that might be expected for 
rapidly diverging flows (Reynolds 1986). For consistency purposes, 
we have used the same archival data as in AMI Consortium: Scaife 
et al. (2012a) for the SED, who find a' = -0.11 ±0.09. However, 
we find a' = 0.70 ±0.06 in this work. We attribute this difference 
to the fact that the data is fit with multiple greybody components 
in AMI Consortium: Scaife et al. (2012a). L1527 has a known disc 
(Loinard et al. 2002; Tobin et al. 2010) which dominates the SED 
at sub-mm frequencies, while the envelope contribution dominates 
at far-infrared frequencies. The single greybody model used in this 
work does not provide the best fit to both the envelope and disc. 
However, in both scenarios the SED model here predicts /3 ^ 1, 
consistent with the presence of large dust grains in the disc (AMI 
Consortium: Scaife et al. 2012a). 

Overall, the values for the spectral indices are consistent as 
all sources listed in Table 5 and the protostellar sources listed in 
Table 8 have known outflows. In the cases of NGC 2264 G and 



L723, higher spatial resolution data can help constrain the emission 
mechanism for the driving sources of these outflows and identify 
non-thermal components. 



6.6 Correlations 

Radio luminosity has been shown to correlate with global prop- 
erties such as bolometric luminosity, envelope mass, and outflow 
force (AMI Consortium: Scaife et al. 2011a,b, 2012b). We do not 
repeat these discussions here, but combine our data with the previ- 
ous work to redraw the correlations, see Figures 3 and 4. We use 
the radio luminosity (Sjad.ie*^ ) results calculated by subtracting the 
predicted greybody contribution at 16 GHz (5gb,i6) from the mea- 
sured flux density (SkjGHz)' in the correlations that follow to ensure 
that the values used are representative of only the outflow compo- 
nent of the radio emission. We compare the correlations between 
both SED scenarios and examine the relation of this sample to ex- 
isting correlations. 



6. 6. 1 Note on Classification 

AMI Consortium: Scaife et al. (2011b) found that when compar- 
ing bolometric temperatures determined using IRAS data (Hatchell 
et al. 2007b) and those found from fluxes calcuated using Spitzer 
SEDs (Dunham et al. 2008), the IRAS temperatures were, on av- 
erage, a factor of 1.6 higher than the corresponding Spitzer val- 
ues. Bolometric temperature is one of the indicators for protostellar 
class and within this sample, one of the sources with 70 ^ Tboi ^ 
150 (HH 111) will change from Class I to when the Spitzer tem- 
perature is used. 



6.6.2 Correlation with bolometric luminosity 

Figure 3 shows the correlation of radio luminosity with bolomet- 
ric luminosity for (a) Scenario (i), where T^ is fixed, and (b) Sce- 
nario (ii), where T^ varies. These new data extend the trend into 
higher luminosities. We note that the bolometric luminosity for 
the Serpens MMS I source used here is Lboj = 17.3 Lq (Enoch 
et al. 2009a), however significantly higher luminosities have been 
reported for this object in the literature which could explain the 
offset from the trend (see §5.9). The value of a' is unconstrained 
for HH 26 IR due to insufficient low-frequency data. As a result, 
our method calculates a significant value for S^bid, resulting in a 
value of 5'rad,l6'^^ consistent with zero due to the low flux density 
detected by AMI. HH 26 IR is therefore not included in the plots. 
Data below 8 GHz is needed to constrain the spectral index of this 
source and establish an accurate measure of the radio luminosity. 

Two sources, L1551 NE and HH 25 MMS, have their loca- 
tion change significantly between Figures 3(a) and (b), due to the 
large difference in their radio luminosities between the two SED 
scenarios (as discussed extensively in §6). Larger radio luminosi- 
ties were calculated for both sources in the case where T^ is fixed. 
However, the lower radio luminosity calculated from the scenario 
where Tj is allowed to vary for HH 25 MMS is in better agreement 
with the trend than the higher value from the scenario where Tj is 
fixed. The opposite is true for LI 551 NE, and the radio luminos- 
ity calculated for this source from Scenario (i) is in better agree- 
ment with the trend than the one calculated from Scenario (ii). Like 
HH 26 IR, L1551 NE deviates from the trend due to a radio lumi- 
nosity Sfadigrf ~ when T^ varies, and more observations at these 
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frequencies are needed to constrain tlie spectral index and establish 
the radio luminosity. 

We use a two sample Kolmogorov-Smimov (KS) test to de- 
termine whether the sample from this paper is consistent with the 
known correlation. We reject the null hypothesis of no difference 
between the datasets if the probability is less than the standard cut- 
off value of p=0.05. We find p-values of p=0.021 and =0.037 for 
Scenarios (i) and (ii) respectively, indicating that the sample from 
this work is statistically different from the previous sample in both 
scenarios. 

We attribute this discrepancy to the values of the bolometric 
luminosity, rather than the radio luminosity, as no similar effect 
is seen in the correlation with envelope mass (see § 6.6.3). There 
are a number of possible explanations. The data set from previous 
papers (open circles in Fig. 3) has bolometric luminosity values 
dominated by a single catalogue (Hatchell et al. 2007b), whereas 
the current sample has values from a much wider range of refer- 
ences. We might therefore expect a larger degree of scatter in the 
current dataset due to differences in methodology. Alternatively, 
it is possible that the bolometric luminosities derived in Hatchell 
et al. (2007b) which dominate the archival dataset have themselves 
a systematic offset from those found more generally in the litera- 
ture, however this seems unlikely as it has not been noted in previ- 
ous comparisons. A further possibility is that the current sample is 
(on average) biased towards sources with high ratios of bolometric 
luminosity to envelope mass. Since these sources are preferentially 
selected to be well known radio protostars it is likely that this effect 
is present at some level. We therefore suggest that the underlying 
cause of the discrepancy is a combination of both this last possibil- 
ity and the first. 
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6.63 Correlation with envelope mass 

Figure 4a shows the correlation between radio luminosity and en- 
velope mass. The envelope masses listed in Column [8] of Table 3 
were determined using the 850 fim flux density from the Funda- 
mental Map Object Catalogue (di Francesco et al. 2008) in the fol- 
lowing way: 



M = 



Syd^ 



KvBviTi) 



(6) 



where Ky is the opacity index, Bv{T^) is the Planck function at a 
dust temperature T^, and the scaling factor [K:v6v(?d)] uses the 
scaled value from Hatchell et al. (2007b), applied in AMI Consor- 
tium: Scaife et al. (2012b). This fixes us on the same scale as Enoch 
et al. (2009a) and allows for comparisons between the datasets. Dis- 
tances used are listed in Column [3] of Table 3. Froebrich (2005) 
find that lower dust temperature of 15 K e.g., as suggested by Motte 
& Andre (2001) and Hatchell et al. (2007b), would lead to envelope 
masses that are a factor of 1 .6 higher, however this should not affect 
our results as we are maintaining consistency within our sample. 
The SCUBA (di Francesco et al. 2008) 850 /im flux density is un- 
resolved for NGC 1333 IRAS 2A and B, and so we use the values 
from Enoch et al. (2009b) for these objects. 

The data from this paper are combined with those from AMI 
Consortium: Scaife et al. (201 la,b) to examine the correlation of ra- 
dio luminosity with envelope mass. It can be seen from Figure 4(a) 
that the data from this work are consistent with those from the lit- 
erature. A best-fitting correlation using the combined data-set is 



log, 



Lk„i (Lm) 



(b) 

Figure 3. (a) Correlation of bolometric luminosity with 1.9 cm radio lumi- 
nosity when Ti is fixed, and (b) correlation of bolometric luminosity with 
1.9 cm radio luminosity when Ti varies. Filled circles represent Class 
sources from this sample, stars are Class I, and unfilled circles are previous 
data. Best fitting correlations from AMI Consortium: Scaife et al. (2011b) 
are shown as solid lines. 



found to be. 



log[Lrad(tnJykpc )] = 



(-I.83±0.25) 

-KO.96±O.41)log[Menv(M0) 



Following AMI Consortium: Scaife et al. (2011b) we also 
compare the data to a power-law with index 1.33, shown as a 
dashed line in Figure 4(a). We note that the fitted correlation is 
consistent with this index at the 1 CT level. 

Class I sources are expected to have smaller envelope masses 
than Class by definition (see §1) but from these data there is no 
evident deviation from the trend based on evolutionary class. In 
fact, the Class I sources in this sample are all found along the -l-l cr 
level line of the correlation. 

From the KS test we find p-values of p=0.374 and p=0.54 for 
Scenarios (i) and (ii) respectively, indicating that the sample from 
this work is consistent with the correlations found from the previ- 
ous sample in both scenarios. . 
Observational bias in the L^d Menv correlation We note that 
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the correlation fitted directly to the data shown in Figure 4 will 
be subject to an observational bias due to the low significance end 
of the data. This form of Eddington bias causes the slope of the 
fitted correlation to be artificially shallow due to the increased data 
scattered above the correlation, relative to that scattered below the 
correlation which will be missed due to the detection limit. Such 
a bias can be corrected statistically, and following the method of 
Herranz et al. (2006) we use a maximum likelihood estimator for 
the slope, ra, to recover a value of ra = 1.01. 

The correlation with bolometric luminosity will also be sub- 
ject to a similar bias. However, we do not correct the fit in this case 
as the detection limit is not the only criterion for detection of a ra- 
dio counterpart in this case, as shown in AMI Consortium: Scaife 
etal. (2011b). 




6.6.4 Correlation with outflow force 

The 1 .9 cm radio luminosity data presented here are consistent with 
the loose correlation with outflow force measured in AMI Consor- 
tium: Scaife et al. (201 la,b), see Figure 4(b). As discussed in AMI 
Consortium: Scaife et al. (201 la) the lack of a well defined correla- 
tion in this case may by due to the errors inherent in measurements 
of outflow force. Class objects are generally considered to have 
more energetic outflows than Class I, but there is no clear evolu- 
tionary division in the correlation presented here. 

Many studies increase their outflow forces by a factor of 10 
to account for optical depth and inclination effects, following Bon- 
temps et al. (1996). From the current sample, only two the values 
for outflow force from the literature have been corrected in this 
way: L1527 and Serpens MMS (see Table 3). As the majority of 
the data included in Figure 4(b) has this correction implemented, 
we apply the same adjustment to the values in Table 3 for the pur- 
poses of the plot. 
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6.7 Variability 

Radio variability might be expected for Class and I sources due to 
episodic accretion rates and mass ejection in the outflows (Kenyon 
& Hartmann 1995; Young & Evans 2005). However, only two 
sources, L1551 IRS 5 and Serpens SMM 1, show strong signs of 
variability based on examination of their spectral energy distribu- 
tions (see Appendix B). This could be due to the fact that there is 
a comparatively larger amount of archival data for these sources, 
and it is likely that a majority of the sources in this sample pos- 
sess variable natures but more observations at multiple epochs are 
needed to confirm this. In the case of NGC 2264 G, evidence for 
variability has been shown (Gomez et al. 1994) for an object in the 
same region as the driving source of this outflow which will affect 
the results of unresolved data. 

L1551. Rodriguez et al. (2003a) have compiled measurements of 
L1551 IRS 5 at 2 cm from 1983 to 1998 to determine the proper 
motion of the multiple sources, and by comparing different epochs 
they found that the southern component underwent a major ejection 
event during the late 1980s and showed considerable elongation in 
the E-W direction. By comparing data at 5 and 15 GHz taken at 
different epochs (see Appendix C) there is clear evidence of vari- 
ability. The flux densities measured at 5 GHz range by a factor of 
~ 1.7 over all epochs, while the flux densities measured at 15 GHz 
range by a factor of ~ 2.2. In addition, the VLA observations of 
Keene & Masson (1990) were measured in September 1988, and 
their flux density at 15 GHz is a factor of ~ 1.7 higher than the 
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Figure 4. (a) Correlation of envelope mass with 1.9 cm radio luminosity, 
where the solid line indicates the best fitting correlation using the com- 
bined data set, see § 6.6.3, the dashed lines indicate a slope of 1.33 and the 
f (J limits on the distribution of data relative to this slope (see AMI Con- 
sortium: Scaife et al. (2011b)); (b) correlation of 1.9 cm radio luminosity 
with outflow force. Filled circles represent Class sources from this sam- 
ple, stars are Class I, and unfilled circles are previous data. The solid line 
indicates the shock ionisation model of Curiel et al. (1989) with a tempera- 
ture of Te = 10^ K. The upper (lower) dashed line indicates the same model 
with a temperature of T^ = 3000 K (10^ K; AMI Consortium: Scaife et al. 
(2011b)). 



flux density measured in this work at a similar frequency. Lim & 
Takakuwa (2006) detect three circumstellar dust disks in this re- 
gion, and therefore luminosity variations might be expected. These 
authors have discussed that the flux of the N jet (see §5 in this work) 
is typically higher than the flux of the S jet, except during the epoch 
of the major ejection event of the southern jet. 
Serpens. Evidence for variability has been shown for Ser- 
pens SMM 1 by AMI Consortium: Scaife et al. (2012b) and refer- 
ences therein. Curiel et al. (1993) suggest that the knotty structure 
in the NW and SE components is the result of a discrete ejection of 
material ("bullets") from the source, which could account for the 
radio flare event in the archival data AMI Consortium: Scaife et al. 
(2012b) that appears to have faded. The observations at 1.9 cm for 
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SMM 1 presented in this worlc were talcen 11 October 2011 and 
the flux density (see Table 5) is a factor of ~ 1.5 higher than that 
presented in AMI Consortium: Scaife et al. (2012b) who measure a 
flux density of 4.736 ±0.237 at 1.8 cm from observations taken be- 
tween December 2010 and January 2011. Our measured flux den- 
sity is within the uncertainty limits of the 15 GHz flux density of 
10± 3 mjy found by Snell & Bally (1986). Enoch et al. (2009a) de- 
tect a high mass disc around this source which should be unstable 
and therefore exhibit luminosity variations. 

HH 1 The Herbig-Haro object HH 1 has been previously detected 
at centimetre wavelengths (e.g. Rodriguez et al. 1990, 2000) how- 
ever, here we do not detect HH 1 at 1.9 cm. We calculate a flux 
density upper limit of 0.34 mJy for this object, inconsistent with the 
value of 0.78 mJy at 2 cm found by Rodriguez et al. (1990). How- 
ever, Pravdo et al. (1985) also failed to detect HH 1 in the radio, 
suggesting that this source is variable. Indeed multi-wavelength ob- 
servations of this object have been discussed by Raga et al. (201 1), 
and HH 1 has been shown to be variable at radio wavelengths. Fu- 
ture observations with the AMI and eMERLIN at the same epoch 
will provide a better picture of the variability of these objects. 



• We examine correlations between the radio luminosity and 
bolometric luminosity, envelope mass, and outflow force and find 
that these data are broadly consistent with correlations found by 
previous samples. We find a slight bias in the correlation with bolo- 
metric luminosity towards higher values of Lboi- 

• Evidence for variability has been shown for L1551 IRS 5 and 
Serpens SMM 1 by close inspection of their SEDs. It is apparent 
from their histories that they have undergone radio flare events, 
supporting non-steady accretion accretion models of proto stellar 
evolution. We also find variability in the Herbig-Haro object HH 1. 
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7 CONCLUSIONS 

We have presented 16 GHz observations made with AMI of a sam- 
ple of low-mass young stars driving jets. 

• The data presented in this work provides support for free-free 
emission as the dominant mechanism for producing the radio lumi- 
nosity from these YSOs. Across the AMI channel bands 80 per cent 
of the detected protostellar sources in the target sample have spec- 
tral indices —0.1 ^ OfAMl ^ 2 consistent with free-free emission 
and we find an average spectral index «ami = 0.42 ± 0.59. Al- 
though the flux densities integrated over the whole of the source 
are in general indicative of free-free emission, it could be possi- 
ble that subcomponents within the source show non-thermal (i.e. 
negative) spectra. 

• We have also presented and examined spectral energy distri- 
butions for each detected protostellar source, combining the AMI 
data with available archival data from an extensive literature search. 
We tested two scenarios to obtain maximum likelihood values for 
the spectral indices, dust opacity indices, and dust temperatures. In 
Scenario (i) we fix T^ for each source based on protostellar evo- 
lutionary class and find that 80 per cent of the objects in the target 
sample have spectral indices consistent with free-free emission, and 
we calculate an average spectral index a' = 0.20 ± 0.41 consistent 
with the value for a well coUimated outflow. In Scenario (ii) we 
do not fix Tj, and we find the same average value for the spectral 
index. These values for the spectral indices are consistent with ex- 
pectations as all target sources have known outflows. 

• We examine the errors associated with determining the radio 
luminosity and find that the dominant source of error is the uncer- 
tainty on j3. 

• In both scenarios we find that NGC 2264 G and L723 
have a' ^ —0.1, inconsistent with free-free emission. For these 
two sources there are other objects within close proximity (i.e. 
within the synthesized beam) that have non-thermal spectral in- 
dices that affect the integrated flux densities. Higher spatial resolu- 
tion data can help constrain the emission mechanism for the driving 
sources of these outflows. We also find a' > 2 in some cases (e.g. 
NGC 1333 IRAS 2A) suggesting that the lower frequency data is 
an extension of the greybody emission and that there is no free-free 
component detectable in these data. 
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APPENDIX A: AMI MAPS 

Maps of each field are shown. AMI data is shown uncorrected for 
the primary beam response. The positions of known Class and 
Class I objects are indicated as crosses (-I-) and stars (*), respec- 
tively. The positions for those found in Perseus (L1448 and HH 7- 
11) are from Hatchell et al. (2007b), HH 1-2 from W. J. Fischer 
et al. (2010), and Serpens from Duarte-Cabral et al. (2010). Crosses 
(x) represent known sources of unknown evolutionary class. We 
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plot the half power point of the primary beam as a solid circle 
(PS 6arcmin at 16 GHz) and the FWHM of the PSF as the filled 
ellipse in the bottom left comer (see Table 1). Contours at 5, 10, 
15, 20 CTrms etc, where CTnns values for each source are listed in Ta- 
ble 1. 



APPENDIX B: SPECTRAL ENERGY DISTRIBUTIONS 

The observed radio spectra over the AMI frequency channels 3-8 
combined with flux densities from the literature. Maximum likeli- 
hood results from Scenario (ii) are overlaid (see Table 8). The list of 
archival data used in the spectral energy distributions can be found 
in Appendix C. Only data v < 3 THz (1 > 100 /xm) were included 
in the fit, but IRAS data v > 3 THz are included in the plots for 
illustration. 



APPENDIX C: SED REFERENCE TABLE 

An extensive literature search was conducted for unresolved, inte- 
grated flux densities to include in the spectral energy distributions. 
It should be noted that Wendker (1995) was a useful reference. 
High resolution data that were highly discrepant due to flux loss or 
data with high uncertainties (e.g. 450 flm data from di Francesco 
et al. 2008) were not included. Where uncertainties were not pro- 
vided, an error of 10 per cent was used in the model fittings and 
this is indicated by a ' . 
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Figure Al. The AMI 16 GHz combined-channel map for each source. 
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Figure A2. The AMI 16 GHz combined-channel map for each source - continued. 
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Figure Bl. The spectral energy distribution for each source along with the model fit. 
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Figure B2. The spectral energy distribution for each source along with the model fit - continued. 
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Figure B3. The spectral energy distribution for each source along with the model fit - continued. 
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Table CI. Reference list for the spectral energy distribution data. 



Source V (GHz) 


Sv (Jy) 


Reference 


L1448 IRS 3 5 


1.15xl0-3t 


Curieletal. (1990) 


8 


1.6x10-31- 


Reipurth et al. (2002) 


14.62 


2.40x10-3 ±0.25 


this work 


15 


1.82xl0-3t 


Curiel et al. (1990) 


15.37 


2.34x10-3 ±0.26 


this work 


16 


2.06x10-3 ±0.12 


AMI Consortium: Scaife et al. (201 lb) 


16.12 


2.41x10-3 ±0.37 


this work 


16.87 


2.26x10-3 ±0.24 


this work 


17.62 


2.62x10-3 ±0.29 


this work 


91 


0.1 ±0.01 


Shirley et al. (2011) 


231 


2.6±0.14 


Barsony etal. (1998) 


231 


3t 


Motte & Andre (2001) 


273 


4.7 ±0.05 


Enoch et al. (2006) 


273 


2.3 ±0.20 


Barsony etal. (1998) 


349 


7.12±0.33 


Shirley etal. (2011) 


353 


10.18+ 


di Francesco et al. (2008) 


353 


16.86+ 


Hatchell et al. (2007b) 


353 


8.37±0.91 


Chandler & Richer (2000) 


353 


14.7±0.62 


Shirley et al. (2000) 


375 


7.8 ±0.64 


Barsony etal. (1998) 


400 


10.9±1.15 


Chandler & Richer (2000) 


666 


56.1 ±12.28 


Chandler & Richer (2000) 


666 


100.3 ±12.40 


Shirley et al. (2000) 


857 


91.5±22.78 


Chandler & Richer (2000) 


857 


45±3 


Barsony etal. (1998) 


3xl03 


112±20.15 


Barsony et al. (1998) 


5x10^ 


32±6.12 


Barsony etal. (1998) 


1.2x10" 


5.75 ±1.2 


Barsony etal. (1998) 


2.5x10'' 


0.69 ±0.15 


Barsony et al. (1998) 


HH7-11 5 


7.3xl0-"t 


Rodriguez et al. (1999) 


8 


1.01xl0-3t 


Rodriguez etal. (1999) 


14.62 


3.43x10-3 ±0.20 


this work 


15.37 


3.10x10-3 ±0.18 


this work 


16.12 


3.91x10-3 ±0.22 


this work 


16.87 


3.43x10-3 ±0.19 


this work 


17.62 


3.83x10-3 ±0.24 


this work 


43 


1.08x10-2 ±0.06 


Anglada et al. (2004) 


100 


7.89 ±0.79 


Chen et al. (2009) 


273 


4.46t 


Enoch et al. (2006) 


353 


IO.3I" 


di Francesco et al. (2008) 


353 


12.7+ 


Hatchell et al. (2007b) 


353 


14.9±1.2 


Chandler & Richer (2000) 


400 


21.5±2.8 


Chandler & Richer (2000) 


666 


119±24 


Chandler & Richer (2000) 


666 


159.9'^ 


Hatchell et al. (2007b) 


857 


203 ±61 


Chandler & Richer (2000) 


3x10^ 


381 ±23 


IRAS 


5x10^ 


204 ±20 


IRAS 


1.2x10" 


46.5 ±2.8 


IRAS 


2.5x10" 


13.6±3.7 


IRAS 


L1551IRS5 1.5 


2.28xl0-3t 


Bieging & Cohen (1985) 


1.5 


2.8x10-3 ±0.9 


Snell et al. (1985) 


1.7 


4xl0-3t 


Rodriguez etal. (1989a) 


5 


4.2x10-3+ 


Duncan et al. (1987) 


5 


3.5x10-3 ±0.5 


Cohen et al. (1982) 


5 


3xlO-3t 


Bieging et al. (1984) 


5 


4.69xl0-3t 


Bieging & Cohen (1985) 


5 


4.3x10-3 ±0.5 


Snell et al. (1985) 


5 


5x10-3 ±0.5 


Evans etal. (1987) 


5 


4.1x10-3 ±0.4 


Keene & Masson (1990) 


8 


4.7x10-3 ±0.5 


Keene & Masson (1990) 


14.62 


4.74x10-3 ±0.25 


this work 
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Source 



V (GHz) Sv (Jy) 



Reference 



L1551IRS5cont. 



L1527 



15 


2.25xl0-3t 


Rodriguez et al. (1986) 


15 


4.59xl0-3t 


Bieging & Cohen (1985) 


15 


4.9x10-3 ±0.5 


Keene & Masson (1990) 


15 


3.6x10^3 _|_ 0.2 


Rodriguez et al. (2003a) 


15 


2.8x10-3 ±0.2 


Rodriguez et al. (2003a) 


15 


3.3x10-3 ±0.2 


Rodriguez et al. (2003a) 


15 


2.7x10-3 ±0.2 


Rodriguez et al. (2003a) 


15.37 


3.10x10-3 ±0.18 


this work 


16.12 


3.91x10-3 ±0.22 


this work 


16.87 


3.43x10-3 ±0.19 


this work 


17.62 


3.83x10-3 ±0.24 


this work 


22.5 


1.7x10-2 ±0.4 


Torrelles et al. (1985) 


22.5 


7.3x10-3 ±1.5 


Keene & Masson (1990) 


23.7 


7x10-3 ±1.4 


Gomez etal. (1993) 


88 


0.09 ±0.03 


Keene & Masson (1990) 


90 


0.12+ 


Altenhoff etal. (1994) 


98 


0.13 ±0.004 


Ohashi etal. (1991) 


110 


0.13 ±0.03 


Keene & Masson (1990) 


112 


0.15+ 


Wilking et al. (1989) 


113 


0.17t 


Keene & Masson (1986) 


230 


1.5t 


Cabrit & Andre (1991) 


230 


1.57±0.02 


Reipurth et al. (1993) 


231 


3.4+ 


Motte & Andre (2001) 


240 


2.37 ±0.48 


Keene & Masson (1990) 


250 


1.72+ 


Altenhoff etal. (1994) 


300 


5.7±1.3 


Keene & Masson (1990) 


345 


6.36 ±0.06 


Reipurth et al. (1993) 


353 


19.01+ 


di Francesco et al. (2008) 


353 


12.1 ±1 


Chandler & Richer (2000) 


400 


18.2±2.4 


Chandler & Richer (2000) 


666 


94 ±19 


Chandler & Richer (2000) 


857 


164 ±49 


Chandler & Richer (2000) 


3x10' 


458+ 


Reipurth et al. (1993) 


5xl03 


373+ 


Reipurth et al. (1993) 


1.2x10"' 


106+ 


Reipurth et al. (1993) 


2.5x10"' 


10+ 


Reipurth et al. (1993) 


5 


6.8x10-** ±0.4 


Melis etal. (2011) 


8.5 


8.10x10-'' ±0.3 


Melis etal. (2011) 


14.62 


1.04x10-3 ±0.11 


this work 


15.37 


1.04x10-3 ±0.14 


this work 


16 


0.9x10-3 ±0.03 


AMI Consortium: Scaife et al. (2012a) 


16.12 


1.12x10-3 ±0.13 


this work 


16.87 


1.20x10-3 ±0.16 


this work 


17.62 


1.38x10-3 ±0.15 


this work 


22.5 


1.4x10-3 ±0.1 


Melis etal. (2011) 


43.3 


4.4x10-3 ±0.6 


Melis etal. (2011) 


HI 


4.7x10-2 ±0.56 


Ohashi et al. (1997) 


230 


1.5"^ 


Motte & Andre (2001) 


230 


0.35+ 


Gramajo et al. (2010) 


353 


11.78+ 


di Francesco et al. (2008) 


353 


0.90 ±0.01 


Gramajo et al. (2010) 


375 


0.50+ 


Gramajo et al. (2010) 


666 


2.85 ±0.22 


Gramajo et al. (2010) 


666 


3.21+ 


Gramajo et al. (2010) 


857 


12+ 


Gramajo et al. (2010) 


1.87x103 


68.8+ 


Gramajo et al. (2010) 


3x10^ 


71.3+ 


Gramajo et al. (2010) 


3x103 


72+ 


Gramajo et al. (2010) 


3x103 


73.3±11.7 


IRAS 


5xl03 


17.4+ 


Gramajo et al. (2010) 


5x10' 


18+ 


Gramajo et al. (2010) 


5x10^ 


17.8±1.6 


IRAS 


1.2xl03 


0.74 ±0.07 


IRAS 
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Table C3. Reference list for the spectral energy distribution data - continued. 



Source V (GHz) 


Sv (Jy) 


Reference 


HH1-2MMS1 1.5 


6xl0-''±2 


Pravdo et al. (1985) 


1.5 


S-exlO-*"! 


Rodriguez etal. (1990) 


5 


1.19x10"' ±1.6 


Morgan etal. (1990) 


5 


1.2x10"' ±0.04 


Pravdo et al. (1985) 


5 


1.17xKr't 


Rodriguez etal. (1990) 


14.62 


1.30x10-' ±0.28 


this work 


15 


1.54x10-' ±0.18 


Pravdo et al. (1985) 


15 


1.75xlO-'t 


Rodriguez etal. (1990) 


15.37 


xUr' ±0.18 


this work 


16.12 


xUr' ±0.22 


this work 


16.87 


xlO-'±0.19 


this work 


230 


0.65 ±0.02 


Reipurth etal. (1993) 


273 


0.92 ±0.03 


Dent etal. (1998) 


345 


1.67 ±0.02 


Reipurth etal. (1993) 


353 


9.44+ 


di Francesco et al. (2008) 


354 


1.2 ±0.24 


W.J. Fischer etal. (2010) 


375 


2.65 ±0.03 


Dent etal. (1998) 


666 


18.9±0.35 


Dent etal. (1998) 


857 


33.95 ±0.29 


Dent etal. (1998) 


857 


± 


W.J. Fischer etal. (2010) 


1.88xl03 


75.7± 15.14 


W.J. Fischer etal. (2010) 


4.29x10' 


26.6 ±2.66 


W.J. Fischer etal. (2010) 


HH 26 IR 8 


3.8x10-* ±0.6 


Anglada etal. (1998) 


8 


1.4x10-'* ±0.3 


Gibb (1999) 


16.12 


3.92x10-'' ±0.76 


this work 


231 


0.32+ 


Lis et al. (1999) 


353 


0.97+ 


di Francesco et al. (2008) 


353 


0.9 ±0.09 


Nutter & Ward-Thompson (2007) 


666 


1.4 ±0.28 


Nutter & Ward-Thompson (2007) 


857 


6.3* 


Lis et al. (1999) 


3x10' 


67.93* 


Antoniucci et al. (2008) 


4.29x10' 


7.9* 


Antoniucci et al. (2008) 


5x10' 


20.87* 


Antoniucci et al. (2008) 


1.2x10^ 


5.13+ 


Antoniucci et al. (2008) 


1.25x10' 


1.8* 


Antoniucci et al. (2008) 


2.5x10'' 


1.93+ 


Antoniucci et al. (2008) 


HHlll 5 


8.3x10-* ±0.9 


Rodriguez et al. (2008b) 


8 


9.5x10-* ±0.4 


Rodriguez & Reipurth (1994) 


14.62 


3.11x10-' ±0.53 


this work 


15 


2.3x10-' ±0.32 


Rodriguez & Reipurth (1994) 


15.37 


2.68x10-' ±0.28 


this work 


16.12 


2.39x10-' ±0.26 


this work 


16.87 


2.86x10-' ±0.30 


this work 


17.62 


2.51x10-' ±0.34 


this work 


43 


5.15x10-' ±0.52 


Rodriguez et al. (2008b) 


110 


4.6x10-2 ± 0.46 


Stapelfeldt & Scoville (1993) 


230 


0.49 ±0.02 


Reipurth etal. (1993) 


230 


0.47 ±0.05 


Wilking et al. (1989) 


273 


0.75 ±0.01 


Dent etal. (1998) 


285 


0.69+ 


Stapelfeldt & Scoville (1993) 


345 


1.39 ±0.03 


Reipurth et al. (1993) 


353 


3.09+ 


di Francesco et al. (2008) 


375 


1.8 ±0.02 


Dent etal. (1998) 


666 


8.22 ±0.26 


Dent etal. (1998) 


857 


14.2±0.9 


Dent etal. (1998) 


3x10' 


71.1+ 


Reipurth etal. (1993) 


5x10' 


43.5+ 


Reipurth etal. (1993) 


1.2x10^* 


6.8* 


Reipurth et al. (1993) 


2.5x10^* 


0.3* 


Reipurth etal. (1993) 


NGC2264G 1.5 


8x10-" ±2 


Rodriguez & Curiel (1989) 


5 


1.9x10-' ±0.2 


Rodriguez & Curiel (1989) 


5 


6.7x10-"+ 


Gomez etal. (1994) 
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Source V (GHz) 


Sv (Jy) 


Reference 


NGC 2264 G cont. 8 


5.5x10-*+ 


Gomez etal. (1994) 


16.12 


2.96x10-"' ±2.16 


this work 


273 


0.25 ±0.04 


Ward- Thompson et al. (1995) 


353 


0.71"!" 


di Francesco et al. (2008) 


375 


0.7 ±0.04 


Ward-Thompson et al. (1995) 


666 


3.9±0.1 


Ward-Thompson et al. (1995) 


857 


7.1 ±0.28 


Ward-Thompson et al. (1995) 


3xl03 


17±3 


Ward-Thompson et al. (1995) 


5xl03 


6±1 


Ward-Thompson et al. (1995) 


Serpens SMM 1 1.5 


4.5x10"^ ±0.5 


Condon et al. (1998) 


1.5 


9.5x10-3 ±0.49 


Curiel et al. (1993) 


5 


9.5x10-3 ±0.8 


Snell& Bally (1986) 


5 


7.9x10-3 ±0.48 


Rodriguez et al. (1989b) 


5 


7.6x10-3 ±0.3 


Curiel et al. (1993) 


5 


2.2xl0-3t 


McMullin et al. (1994) 


8 


7.5x10-3 ±0.17 


Curiel et al. (1993) 


8 


7.54xl0-3t 


Eiroa et al. (2005) 


14.62 


8.21x10-3 ±0.77 


this work 


15 


1x10-2 ±0.3 


Snell& Bally (1986) 


15 


8.3x10-3 ±0.3 


Curiel et al. (1993) 


15 


6.2x10-3 ±0.0.44 


Rodriguez et al. (1989b) 


15 


4.8xlO-3t 


McMullin et al. (1994) 


15.37 


7.66x10-3 ±0.0.55 


this work 


16 


4.74 ±0.24 


AMI Consortium: Scaife et al. (2012b) 


16.12 


7.04x10-3 ±0.43 


this work 


16.87 


6.77x10-3 ±0.44 


this work 


17.62 


7.30x10-3 ±0.54 


this work 


23 


5.3x10-3+ 


McMullin et al. (1994) 


43 


1.42x10-2 ±0.04 


Choi (2009) 


88 


0.2+ 


Hogerheijde et al. (1999) 


94 


0.2+ 


Hogerheijde et al. (1999) 


97 


0.14+ 


McMullin et al. (1994) 


111 


0.41+ 


Hogerheijde et al. (1999) 


214 


2.65+ 


Hogerheijde et al. (1999) 


240 


2.3+ 


McMullin et al. (1994) 


273 


3.47 ±0.1 


CasaU etal. (1993) 


353 


15.23t 


di Francesco et al. (2008) 


353 


6.1+ 


Davis et al. (1999) 


666 


35.7+ 


Davis etal. (1999) 


1.88x10^ 


430+ 


McMullin et al. (1994) 


3x10^ 


380+ 


McMullin et al. (1994) 


5x10^ 


152.92+ 


McMullin et al. (1994) 


6x10^ 


88.6+ 


McMullin et al. (1994) 


1.2x10"' 


3.17+ 


McMullin et al. (1994) 


1.5x10" 


2.6+ 


McMullin et al. (1994) 


2.5x10"' 


0.25+ 


McMullin et al. (1994) 


L723 5 


7.5x10-"+ 


Anglada et al. (1996) 


8 


7.6x10-"+ 


Anglada et al. (1996) 


14.62 


5.82x10-" ±0.86 


this work 


15.37 


5.16x10-" ±0.64 


this work 


16.12 


6.53x10-" ±0.78 


this work 


16.87 


6.55x10-" ±0.69 


this work 


17.62 


5.22x10-" ±1.00 


this work 


230 


0.37+ 


Motte & Andre (2001) 


230 


0.36 ±0.02 


Reipurth et al. (1993) 


300 


1±0.5 


Reipurth et al. (1993) 


353 


1.79±0.11 


Shiriey et al. (2000) 


353 


2.24+ 


di Francesco et al. (2008) 


666 


8.5±2.1 


Shiriey et al. (2000) 


750 


13±3 


Reipurth et al. (1993) 


857 


11.3±1.8 


Wu et al. (2007) 


1.54x10^ 


35±7 


Reipurth et al. (1993) 


1.81x10^ 


40 ±12 


Reipurth et al. (1993) 
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Table C5. Reference list for the spectral energy distribution data - continued. 



Source V (GHz) 


Sv (Jy) 


Reference 


L723 cent. 2.08x10^ 


33 ±10 


Reipurthetal. (1993) 


2.14x10' 


23±8 


Reipurthetal. (1993) 


2.31x10' 


32±11 


Reipurthetal. (1993) 


3x10' 


20.7 ±1.7 


Shirley et al. (2000) 


3.16x10' 


27±6 


Shirley et al. (2000) 


5x10' 


6.93 ±0.62 


Shirley et al. (2000) 


1.2x10* 


0.38 ±0.05 


Shirley et al. (2000) 


2.5x10" 


0.28 ±0.06 


Shirley et al. (2000) 


L1251A 5 


1.7x10"* ±0.03 


Beltran et al. (2001) 


8 


2.9x10"*+ 


Reipurth et al. (2004b) 


8 


4.7x10-* ±0.03 


Beltran et al. (2001) 


14.62 


8.57x10-* ±1.77 


this work 


15.37 


9.85x10-* ±1.36 


this work 


16.12 


1.00x10-' ±0.20 


this work 


16.87 


7.61x10-* ±1.58 


this work 


17.62 


1.10x10-' ±0.23 


this work 


230 


0.23 ±0.02 


Rosvick & Davidge (1995) 


273 


0.38 ±0.03 


Rosvick & Davidge (1995) 


353 


5.68"^ 


di Francesco et al. (2008) 


375 


0.71±0.11 


Rosvick & Davidge (1995) 


857 


20.8 ±3.2 


Wu et al. (2007) 


3x10' 


78.5 ±12.6 


Rosvick & Davidge (1995) 


3x10' 


79+ 


Sato & Fukui (1989) 


5x10' 


67.7 ±6.1 


Rosvick & Davidge (1995) 


5x10' 


66+ 


Sato & Fukui (1989) 


1.2x10* 


28.3±1.4 


Rosvick & Davidge (1995) 


1.2x10* 


26t 


Sato & Fukui (1989) 


2.5x10* 


6.2±0.3 


Rosvick & Davidge (1995) 


2.5x10* 


5' 


Sato & Fukui (1989) 


L1448 C 8 


2.3x10-*+ 


Reipurth et al. (2002) 


15 


5.6x10"* ±0.5 


Curiel et al. (1990) 


16 


5.37x10-* ±0.33 


AMI Consortium: Scaife et al. (201 lb) 


16.12 


7.58x10-* ±0.94 


this work 


86 


2.6x10-2 ±0.2 


Shirley et al. (2000) 


87 


1.2x10-2+ 


Bachiller et al. (1991) 


87 


1.6x10-2+ 


Guilloteauetal. (1992) 


91 


3.92x10-2 ±0.39 


Shirley etal. (2011) 


115 


9.1xl0-2±0.2 


Shirley et al. (2000) 


230 


1±0.1 


Barsony etal. (1998) 


230 


0.74±0.11± 


Shirley et al. (2000) 


230 


0.91+± 


Motte & Andre (2001) 


230 


0.58+ 


Bachiller et al. (1991) 


273 


2.04+ 


Enoch et al. (2006) 


273 


1±0.1 


Barsony etal. (1998) 


349 


1.97 ±0.13 


Shirley etal. (2011) 


353 


2.98+m 


di Francesco et al. (2008) 


353 


6.51 + 


Hatchell et al. (2007b) 


353 


3.95 ±0.24 


Shirley et al. (2000) 


353 


5.34 ±0.43 


Chandler & Richer (2000) 


375 


3 ±0.3 


Barsony etal. (1998) 


400 


6.91 ±0.91 


Chandler & Richer (2000) 


666 


68.4+ 


Hatchell et al. (2007b) 


666 


21±2 


Barsony etal. (1998) 


666 


31.8±5.5 


Shirley et al. (2000) 


666 


34.2 ±6.9 


Chandler & Richer (2000) 


857 


58±18 


Chandler & Richer (2000) 


857 


30±3 


Barsony etal. (1998) 


3x10' 


70.3 ±14.8 


Barsony et al. (1998) 


5x10' 


31.2±6.5 


Barsony etal. (1998) 


1.2x10* 


2.9 ±0.6 


Barsony etal. (1998) 


2.5x10* 


0.33 ±0.07 


Barsony etal. (1998) 


NGC 1333 IRAS 2A 5 


6x10-5 ±1 


Rodriguez et al. (1999) 
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Source V (GHz) 


Sv (Jy) 


Reference 


NGC 1333 IRAS 2A cont. 8 


2.2xl0-''1" 


Reipurth et al. (2002) 


8 


2.5x10-'' ±0.4 


Rodriguez etal. (1999) 


16 


3.2xl0-''±0.25 


AMI Consortium: Scaife et al. (2011b) 


16.12 


3.92x10-'' ±0.96 


this work 


43 


lxl0-2±0.5 


Anglada et al. (2004) 


86 


3.5xlO-2t 


J0rgensen et al. (2004) 


89 


4x10-2+ 


J0rgensen et al. (2004) 


111 


8.28x10-2 ±0.4 


Looney et al. (2000) 


150 


0.32 ±0.07 


Sandell et al. (1994) 


230 


0.88t 


Motte & Andre (2001) 


230 


0.88±0.12 


Sandell et al. (1994) 


273 


1.46±0.12 


Sandell et al. (1994) 


353 


7.78+ 


Kirk et al. (2007) 


353 


4.79 ±0.39 


Chandler & Richer (2000) 


375 


4.08 ±0.05 


Sandell et al. (1994) 


375 


3.75 ±0.15 


Sandell et al. (1994) 


400 


5.54±0.21 


Sandell et al. (1994) 


400 


6.61 ±0.86 


Chandler & Richer (2000) 


666 


43±11 


Chandler & Richer (2000) 


666 


23.7 ±0.9 


Sandell et al. (1994) 


857 


37.2±1.2 


Sandell et al. (1994) 


857 


74 ±22 


Chandler & Richer (2000) 


3x10^ 


300+ 


Jennings etal. (1987) 


6x10^ 


104+ 


Jennings etal. (1987) 


NGC 1333 IRAS 2B 5 


2.3x10-"' ±0.1 


Rodriguez etal. (1999) 


8 


4x10-" ±0.2 


Rodriguez etal. (1999) 


8 


3.7x10-" ±0.37 


Reipurth et al. (2002) 


16 


1.21x10-3 ±0.05 


AMI Consortium: Scaife et al. (2011b) 


16.12 


5.05x10-" ±1.01 


this work 


43 


5.2x10-3 ±0.3 


Anglada et al. (2004) 


86 


1.2x10-2+ 


J0rgensen et al. (2004) 


89 


1.4x10-2+ 


J0rgensen et al. (2004) 


111 


2.77x10-2 ±0.32 


Looney et al. (2000) 


353 


1.19±0.1 


Chandler & Richer (2000) 


400 


1.53±0.21 


Chandler & Richer (2000) 


666 


7.7 ±1.6 


Chandler & Richer (2000) 


857 


13.4±1.3 


Chandler & Richer (2000) 


L1551 NE 8 


6.6x10-"+ 


Reipurth et al. (2002) 


16.12 


7.43x10-" ±0.92 


this work 


230 


1.5^ 


Motte & Andre (2001) 


230 


0.85 ±0.01 


Andrews & Wilhams (2005) 


230 


0.85 ±0.08 


Moriarty-Schieven et al. (2000) 


238 


0.8±0.12 


Barsony & Chandler (1993) 


272 


1.07±0.11 


Barsony & Chandler (1993) 


353 


2.78+ 


Moriarty-Schieven et al. (2006) 


353 


10.5+ 


di Francesco et al. (2008) 


379 


2.22 ±0.37 


Barsony & Chandler (1993) 


666 


8.33+ 


Moriarty-Schieven et al. (2006) 


677 


16.2±4.6 


Barsony & Chandler (1993) 


857 


22.83 ±0.72 


Andrews & Wilhams (2005) 


3x10^ 


130+ 


Emerson etal. (1984) 


5x10' 


80' 


Emerson etal. (1984) 


1.2x10" 


16+ 


Emerson et al. (1984) 


2.5x10'' 


1.2+ 


Emerson etal. (1984) 


HH 25 MMS 5 


1.5x10-" ±0.4 


Bontemps etal. (1995) 


8 


2.5x10-" ±0.6 


Bontemps et al. (1995) 


8 


1.6x10-" ±0.3 


Gibb(1999) 


16.12 


1.23x10-3 ±0.10 


this work 


231 


1+ 


Lis et al. (1999) 


353 


9.11+ 


di Francesco et al. (2008) 


353 


4+ 


Nutter & Ward-Thompson (2007) 


666 


10.8+ 


Nutter & Ward-Thompson (2007) 
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Table C7. Reference list for the spectral energy distribution data - continued. 



Source 


V (GHz) 


Sv (Jy) 


Reference 


HH 25 MMS cont. 


666 


8.7^ 


Phillips et al. (2001) 




353 


1.4^ 


Phillips et al. (2001) 




857 


24^ 


Lis etal. (1999) 




3x10^ 


70.66"^ 


IRAS 




5x1 ()3 


28.04^ 


IRAS 




5x1 ()3 


32' 


IRAS 




1.2x10* 


1.39+ 


Gezari et al. (1999) 




2.5x10* 


0.51+ 


Gezari et al. (1999) 



